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ABSTRACT 


A procedure for measuring the viscous damping of rela- 
tively large plate material (up to 40 inches x 14 inches x 2 
inches) was developed utilizing the Hewlett-Packard 5451C 
Fourier Analyzer and impulse hammer technique under very low 
stress conditions. Testing environment can be lab air or 
nondistilled water in the temperature range from 30° F to 
Bio? FE. 

The test procedure includes modal analysis that is 
expandable to other geometric shapes and varied material 
such as high damping alloys and composites both metallic 
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fe Ne RODUCEION 


A. GENERAL 

Ship silencing continues to be a major design requirement 
in the construction and operation of submarines and surface 
ships. It has been long established that both equipment and 
personnel are adversely affected by unwatned noise and vibra- 
tion. In today's environment both subsurface and surface, 
with recent advances in acoustic devices the very surviva- 
bility of the platform is directly related to own ship's 
noise. Weapon platforms must be quiet enough to escape detec- 
tion by sophisticated passive sonar devices as well as quiet 
enough to prevent own ship's noise from interfering with 
detection and prosecution of enemy targets. 

The more traditional approach to ship silencing has been, 
and continues to be, that of vibrations isolation. This ap- 
proach requires that all sources of vibration to be placed 
on resilient energy-absorbing mounts. This results ina 
Significant reduction in the transmission of vibrational 
energy from the equipment and main engines to the surrounding 
environment but is not totally satisfactory because these same 
mounts must serve as shock mounts. In addition, any resilient 
mount will have a peak efficiency in particular frequency 
range with decreasing efficiency both above and below that 


frequency range. The optimal solution would be to use 


eZ 





resilient mounts in conjunction with other methods of 
vibrational control to reduce the amount of vibration and 
noise generated as well as reducing the amount transmitted. 

As with any design, today's ships are a compromise of 
requirements and cost considerations. However, the trend 
recently has been toward the design of equipment that is 
inherently quieter, which results in less noise and vibra- 
tional energy at the source. This current trend in design 
is currently being approached from many directions. Among 
these considerations are: 

1. Extreme care in balancing of rotating machinery with 
consideration of vibrations during start up and coast 
down as well as steady state operation. 

2. Selection of operating frequency as far removed from 
resonant frequencies as possible. 

3. Close attention to component tolerances. 

4. Possible use of high damping materials for machine 
elements as well as machine caSings and load bearing 
SEemuceunes:. 

The last category, use of high damping materials, is the 
furthest from being utilized to its maximum potential for 
noise and vibrational reductions. 

Only recently has the designer even considered the inter- 
nal capacity of a material alone with the traditional material 
properties of strength, fatigue resistance, toughness and 


Gorrosion resistance. 
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There are several commercially available high damping 
alloys that meet the strength requirements for most shipboard 
applications. It is recognized that each of the poniidate 
materials may have unique properties that may cast doubt on 
their actual usefulness in the United States Navy. A more 
important problem exists, however. That is, the lack of a 
consistent method of measuring damping for a material under 


low stress, and high frequency range. 


B. DAMPING 

"Structural damping" refers to a structure's or struc- 
tural component's capacity for dissipating energy, or, more 
precisely, to its capacity for removing from a structural 
vibration some of the energy associated with that vibration. 
The energy removed may be converted directly into heat, trans- 
ferred to connected structures or ambient media [Ref. 1]. 

Damping has two primary effects: (1) It limits the 
steady-state motions of structures or systems in situations 
where these motions are controlled by an energy balance; and 
(2) It increases the rates at which the free (i.e., unforced) 
vibrations of structures decay. 

Consider, for example, the classical lumped-parameter 
mMass-spring dashpot system driven by a steady sinusoidal force 
that acts on the mass. For such a system it is possible to 
make the following observations: (1) For excitation frequen- 


cies that are considerably lower than the system resonance 
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frequency the applied force essentially is balanced by the 
Spring force. The mass and dashpot here have virtually no 
effect. (2) For excitation frequencies that are consider- 
ably higher than the system resonance frequency, the applied 
force essentially is opposed by only the inertia force of the 
mass. In this case the spring and dashpot have virtually no 
effect. 

(3) For excitation frequency near the system resonance 
frequency, the spring force and the inertia force essentially 
cancel each other, leaving only the dashpot (i.e., the damping 
element) to oppose the externally applied force. One may 
note that changes in the damping do not affect the system 
response for the first two of the above cases, but for an 
oscillatory force acting at the system's resonance, the vibra- 
tion amplitude decreases with increasing damping. 

In addition to the case of resonant or near-resonant 
excitation, there occur several other important steady-state 
or near-steady-state situations in which the responses of 
systems or structures are controlled by a balance between the 
energy input and the energy dissipations. These situations 
include cases of (1) broad-band excitation, and (2) spatially 
periodic excitation where the spatial period matches that of 
freely propagating waves. Because broad-band excitation gener- 
ally encompasses several structural resonance frequencies, 
the structural responses of the excited modes, and each 


resonant modal response is controlled by damping. 
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If a structure (or a mechanical system) is deflected from 
its equilibrium position and then released, the structure 
vibrates with ever-decreasing amplitude on the Sa aee 
damping, i.e., as the result of energy being removed from 
the oscillatory motions. Greater damping corresponds to the 
dissipation per cycle of a greater fraction of the vibratory 
energy, thus resulting in more tapid decay of the vibrations. 
In a somewhat similar manner, increased damping also results 
in the more rapid decay of freely traveling waves. For exam- 
ple, if a long beam is subjected to an oscillatory transverse 
force of a constant amplitude and frequency, such as near the 
center of the beam, then flexural waves travel away from the 
driving point in both directions. As the result of damping, 
the amplitude of these waves decreases with increasing distance 
from the point of application--with greater damping leading 
to lesser amplitudes at a given distance. 

The practical consequences of effects mentioned above 
generally are the reasons one is interested in systems or 
devices that increase damping. Reductions in resonant or 
random responses result in decreased oscillatory stresses 
and increase in the fatigue life of structures, in the relia- 
bility of mechanical devices, and in mechanical impedance 
(which tends to improve the effectiveness of vibration 
isolation). 

Reduction of the spatially resonant responses of a wall 


Or panel leads to decreased sound transmission through that 


life: 





structure for frequencies above the coincidence frequency. 
Increased attenuation of propagating waves results in lesser 
transmissions of vibrations to neighboring structures. More 
rapid decay of free vibrations reduces the "ringing" sound of 
structures, thus leading to less noise, particularly from 
structures excited by repetitive impacts. This also tends 


to reduce structural fatigue. 


C. MEASURES OF DAMPING 

Damping may be quantified in terms of any of the previously 
discussed primary effects. The corresponding commonly used 
measures of damping, defined below, are interrelated as 


follows: (Refs. 2,3] 





The loss factor n and damping capacity wy are defined 
directly in terms of the cyclic energy dissipation; the 
damping capacity represents the fraction of the system's 
vibrational energy that is dissipated per cycle of the vibra- 
tion, and the loss factor similarly is defined as the frac- 
tion of the system's energy that is dissipated per radian of 
mene vibratory motion. 
are related to the rate of 


On the other hand, T A 


O10)! “aes 
decay of free vibrations. TE denotes the reverberation 


time (in seconds), defined (in analogy to the relatec 


ay, 





room~acoustic measure) as the time within which the vibra- 
tion level of a system vibrating freely at a frequency 

fn (Hz) decreases by 60 dB (i.e., the amplitude decreases to 
1/1000 of its original value). A related measure, decay rate 
A, (dB/sec), represents the rate of reduction of the vibra- 
tion (acceleration or displacement velocity) level. The 
logarithmic decrement 6 is defined as the natural logarithm 
of the ratio of a peak excursion of a freely vibrating system 
to the peak excursion one cycle (period) later. The spatial 
decay A. (dB/wavelength) represents the reduction in the 
steady-state vibration level with distance that occurs along 
a long beam vibrating in flexure. 

It should be noted that none of the measurements of damping 
depend on how the energy is disSipated. Within a cycle these 
measures make no reference to any damping mechanism. On the 
other hand, some other commonly employed measures of damping 
are defined on the basis of viscous damping (1i.e., damping 
that results from a retarding force that is proportional to 
the velocity. The ratio of the magnitude of that force to 
the velocity is called the viscous damping coefficient and 
is commonly designated by c. 

If a simple mass-spring-dashpot system (where the dashpot 
provides a viscous retarding force characterized by c) is 
deflected from equilibrium and released, it typically oscil- 
lates with ever decreasing amplitude. However, if c is made 


large enough, no oscillations occur. Instead, the system 
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creeps toward its equilibrium position, never traversing it. 
The smallest viscous damping coefficient for which this non- 
oscillatory behavior is obtained--i.e., the viscous damping 
that represents the dividing line between oscillatory and 
nonoscillatory behavior is called the critical (viscous) 
damping coefficient ¢c.. It obeys c, = 2vkm where m denotes 
the mass and k the spring stiffness. The "damping ratio" 
C/Cu, also called the fraction of critical damping and often 
given in terms of "percent of critical damping", is widely 
used to indicate damping magnitudes. 

Two other measures of damping are derived from the steady 
state behavior of an ideal linear mass-spring-dashpot system 
that is driven by a Sinusoidal force of constant amplitude 
[Ref. 4]. The amplification at resonance, often called "the 
QO" of the system is defined as the ratio of the amplitude that 
results at resonance to the amplitude that is obtained if 
the force acts quasi-statically (i.e., a frequency consider- 
ably below resonance). The proportional bandwidth b takes 
account of the damping related broadening of the peak ina 
plot of response amplitude versus frequency; this nondimensional 
bandwidth is defined as f/f. where £ denotes the difference 
between the two frequences (one above, and one below the 
resonance frequency f.? at which the square of the response 
amplitude is one-half of its maximum value. For values of 
damping below critical, the aforementioned measure of damping 
are related to each other and to the previously discussed 


loss factor as: 


sis, 





D. DAMPING MECHANISMS 

For linear, viscously damped systems, all of the measures 
of damping discussed previously are independent of the ampli- 
tude. Amplitude independence also occurs for other damping 
mechanisms, and approximate amplitude-independence occurs for 
almost all systems, provided the damping is relatively small. 
Thus, amplitude independence cannot be taken as an indication 
that a system is viscously damped; nevertheless, damping of 
systems with unknown energy dissipation mechanism is often 
characterized in terms of an equivalent viscous damping 
coefficient. 

Indeed, much analysis is carried out with the (usually 
tacit) assumption that the damping is viscous--largely be- 
cause viscous damping leads to linear differential equations 
that can be solved relatively easily. It is fortunate that 
for many practical problems, e.g., where only certain response 
maxima are of concern, the details of the damping mechanism 
are unimportant. Thus, one may obtain reasonable response 
predictions even with inaccurate damping force-versus-velocity 
representations. However, realistic damping models are re- 
quired for the analysis of cases where one is interested in 


details of the response motions. 
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Unlike mass, a single physical phenomenon and stiffness, 
which results from a very few physical effects, damping may 
be caused by a great variety of phenomena. These phenomena 
include mechanical hysteresis (also called material damping 
Or internal friction), electromagnetic effects (notably 
eddy currents), friction due to motion relative to fluids 
or solid surfaces, and energy transport to adjacent structural 
components or fluids (including by acoustic radiation). 

This great variety of phenomena that can produce damping 
generally make damping difficult to predict and to eliminate, 
but enables one to conceive a variety of means for increasing 


ante 
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Il. NATURE OF THE PROBLEM 


A. BACKGROUND 

To date, the majority of damping research has been con- 
ducted on test specimens that are subjected to high stresses 
(1.e., torsional vibrations and vibrating cantilevered beams). 
For most United States Naval applications, structures and 
components are designed for operation at relatively low 
SELesses. 

Damping characteristics are shown to be highly stress 
dependent [Ref. 5]. Therefore, the damping values obtained 
at high stresses do not generally apply for most naval 
applications of ship silencing problems. 

David W. Taylor Naval Ship Research and Development Center 
(NSRDC), Annapolis, Maryland is the focal point for testing 
and evaluating new high damping material for possible use 
for the United States Navy. Realizing the wide scope of the 
work necessary for the testing and evaluation of any new 
material, the author decided to approach NSRDC with a proposal 
to offer assistance in determining the damping of candidate 


material. 


bee UBJECTIVE 
To design a test procedure that allows the measurement of 
damping in a plate specimen at low stress with the following 


variables: 
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l. Plate specimen (40 inches x 14 inches x1 inch) and 
40 inches x 14 inches x 2 inches), or in proportionally 
reduced sizes. 

2. Frequency: Damping testing will be conducted in the 
acoustic range (100-20,000 Hz). 

3. Environment: Testing will be conducted in lab air 
and in nondistilled water. Temperature will range 
Bae bf to JO °F . 

The test procedure includes a complete modal analysis 
that is expandable to other geometric shapes. Actual testing 
and verification of the procedure are conducted on the 
Hewlett-Packard 5451C Digital Fourier Analysis System. 
General information on the HP-5451C system is included in 


Appendix A. 


See oCOrPE OF WORK TO BE COMPLETED 

NSRDC has provided four specimens for testing and to 
establish baseline data for future work. The test specimens 
currently on hand are: A) Cast manganese bronze, code DEQ; 
B) Cast nickel-aluminum bronze, code FTC; C) Steel plate, 
HY-130, code FTW; and D) Aluminum alloy, 5086-H116, code ESX. 
It is anticipated that after baseline data has been compiled 
for both air and water tests that additional specimens will 
be provided by NSRDC. These additional specimens are ex- 
pected to include high damping alloys such as Sonoston and 
Incramute, as well as constrained layer specimens and com- 


posites, both metallic and non-metallic. 
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For the purpose of this thesis the following work was 
completed: (A) design and construction of the test chamber, 
(B) design of a test procedure for the evaluation of damping 
of specimen provided by NSRDC, (C) complete tank character- 
ization by impulse hammer techniques, (D) measurement of 
damping of the cast nickel-aluminum bronze, code FTC, speci- 
men by impulse hammer technique. All of the above testing 


was completed in lab air at normal ambient temperature. 
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Zit. DESIGN OF TEST CHAMBER AND THEORY OF CHARACTERIZATION 


To accomplish the aforementioned objectives four major 
steps must be accomplished: (A) design of a test chamber, 
(B) characterization of the test chamber, (C) design of 
experimental procedure of the measurement of damping, and 
(D) utilization of the above test procedure for damping 


measurement on a plate specimen supplied by NSRDC. 


eee DESoLGN OF TEST CHAMBER 

Design of the test chamber was based on three major 
considerations; (1) size of the largest specimen, (2) weight 
of the largest specimen, and (3) environmental conditions 
desired for testing. In addition, the fixture to support 
the specimen had to hold the specimen rigid so that there 
would be no swinging. 

The overall outside dimension of the chamber (Figure l) 
was chosen to be a rectangular box 48 inches wide by 36 inches 
deep by 72 inches high. This allowed for ample clearance for 
the largest specimen (40 inches x 14 inches x 2 inches) on all 
surfaces and to limit the amount of reflected acoustic energy 
if acoustic absorption material was required on the inside 
surfaces. The test chamber was constructed from one quarter 
inch plate steel supported by a frame constructed from 2 inch 


steel thick wall square tubing. The 2 inch square tubing 
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was welded at all mating edges to give the desired rigidity 
and to prevent water entrapment (future corrosion problems). 

The specimen fixture (Figure 2) was constructed Ree 
3 inches x 2 inches thick wall steel tubing and is designed 
to give a deflection, from a 400 pound specimen, of less than 
OP04 aneh. 

To accommodate various sized specimen the random frequency 
exciter base (Figure 3) can be adjusted within the 2 inch 
rectangular tubing frame work of the upper half of the rear 
panel (Figure 4). The random acoustic exciter can be adjusted 
from 8 inches below the specimen fixture to 24 inches below 
the specimen fixture. This feature allows for a wide variety 
of centered excitation from the end of a specimen plate to 
the center of a specimen plate. 

To allow the removal and placement of various size speci- 
men, the front panel of the test chamber is bolted to the 
test chamber frame and a water gasket is formed by using a 
RTV silicone adhesive sealant. Because of the bonding of 
the RTV all water test must be conducted with the steel front 
panel. For air testing the optional clear front panel may 
be used because no sealant is required. 

The preservation of the test chamber was accomplished by 
sandblasting all steel parts, after protecting front panel 
and specimen fixture bolt threads, one spray paint coat of 
Zinc chromate primer and one hand brushed top coat of enamel. 

To allow testing in lab air and nondistilled water, the 


test chamber is water tight with a 1 inch, valved, drain 
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Inside rear panel of test chamber 
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line near the bottom of the tank. Heating and cooling of the 
Palmers saccomolished by copper tubing (3/8 inch O.D.) coiled 
over entire outside surface of the back panel iEuig anes). 
The coils are spaced 3 inches apart and are 40 inches wide. 
To heat or cool the test chamber a mixture of 50% ethol glycol 
(automative) antifreeze is supplied from a storage tank 
(Figure 6) at the desired temperature and is pumped through 
the heating and cooling coil. 

The entire test chamber is insulated with one and one 
half inches of rigid foam insulation, with the exception of 
the rear panel, which 1s insulated with standard residential 


batt ensulation. 


Pee oenmOny OF THE CHARACTERIZATION OF TEST CHAMBER 

The characterization of the test chamber was accomplished 
by use of the impulse technique for structural frequency 
response testing [Ref. 7] and the Fourier transfer function 
capability of the HP-5451C Fourier Analyzer. 

1. Theory of Frequency Response Function 

The measurement of the frequency response function is 

the heart of modal analysis. The frequency response function 
H(f) is defined in terms of the single input/single output 
system, as the ratio of the Fourier transforms of the system 
output or response v(t) to the system input or excitation 


meee ngquation (3). 


H(f) = (3) 
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Cooling/heating water tank with pump, heater and controller 


Figure 6. 





where: 


V(f) 


II 


Fourier transform of system output v(t) 


iGED) 


Fourier transform of system input u(t). 


The only requirements for a complete description of 
the frequency response function are that the input and output 
Signals be Fourier transformable, a condition that is met by 
all physically realizable systems, and that the input signal 
be non-zero at all frequencies of interest. If the system 
is nonlinear or time-variant, the frequency response function 
will not be unique, but will be a function of the amplitude 
of the input signal in the case of a nonlinear system and 
a function of time in the case of a system with time-varying 
properties. 

The frequency fesponse function may be computed directly 
from the definition as the ratio of the Fourier transforms 
of the output and input signals. However, better results are 
obtained in practice by computing the frequency response 
function as the ratio of the cross-spectrum between the input 
and output to the power spectrum of the input, Equation (4). 
This relationship is derived by multiplying the numerator and 
denominator of the right-hand side of Equation (1) by the 


complex conjugate of the input Fourier transform. 
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where: 


G (f) = U*(f)V(f£), cross-spectrum between u(t) 
uv 7 
and v(t) 
Guy ‘f) = U*(f£f)U(f), power spectrum of u(t) 
U* = complex conjugate of U(f). 


The usefulness of this form of the frequency response function 
can be seen by considering the practical single input/single 
output measurement situation, where m(t) and n(t) represent 
noise at the input and output measurement points, respectively. 

The measured frequency response function H'(f) is 
given by the equation: 


MiGs VT ee ENCE) 
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where the upper case letters denote the Fourier transform of 
the corresponding time domain signals. 

In this form, the measured frequency response will be 
a good approximation of the true frequency response only if 
the measurement noise at both the input and output measure- 
ment points is small relative to the input and output 
Signals. Multiplying the numerator and denominator of the 
right-hand side of Equation (5) by the complex conjugate of 


X(f) yields 


G y ‘f) ~ Gon ‘*) - Gay ‘#) - Gin‘) 
GC, f) + Gm‘? a Gan ce) Gf) 
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Now, if the measurement noise signals m(t) and n(t) 
fee moncoherent with cach other and with the input signal 
u(t), then the expected value of the cross-spectrum terms 


involving m and n in Equation (6) will equal zero, yielding 


Guy ( f) 
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where H(f) is the desired true frequency response function. 

Thus, if the noise-to-signal ratio at the input 
measurement point [G_(£) 7G (£)] is much less than 1, the 
measured frequency response will closely approximate the 
desired true frequency response function. 

It should be pointed out here that there is an inher- 
ent bias error associated with the computation of the cross- 
spectrum and the magnitude of this bias error is inversely 
proportional to the number of averages in the computation. 
Thus, the greater the measurement noise, the greater the num- 
ber of averages required to approach the expected value of 
the cross-spectrum between the input and the output measure- 
ment signals. With measurement techniques employing many 
averages, the bias error can usually be reduced to an insigni- 
ficant level so that it is only necessary to minimize the noise 
in the measurement of the input signal. However, if there is 
Significant measurement noise and only a few averages are 


used, then the computed values of the cross-spectrum terms 
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involving the noise signals in Equation (6) can be large 
relative to the true cross-spectrum, with resulting large 
errors in the measured frequency response function. In 
general, only a few averages are used in the impulse technique; 
otherwise, one of its major advantages--its speed--is lost. 
Therefore, it is important to minimize measurement noise in 
both the input and output Signals when using the impulse 
technique. The cross-spectrum bias error and its effects 

are discussed in more detail in Reference 8. 

Coherence Function. There is another important reason for 
computing the frequency response function in terms cf the 
cross-spectrum: it allows the computation of the coherence 
function between the input and output signals. The coherence 
function is defined by the equation 


Z 
|. GumeCe) | 
(f) = as (8) 


Y xy 6, (£) 6, (4) 

According to the definitions of the power spectrum 
and the cross-spectrum, the coherence function will be iden- 
tically equal to 1 if there is no measurement noise and the 
system is linear. The minimum value of the coherence function, 
which occurs when the two signals are totally uncorrelated, is 
O. Thus, the coherence function 1S a measure of the con- 
tamination of the two signals in terms of noise and nonlinear 
effects, with very low contamination indicated for values 


euose to l. 
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Since the cross-spectrum is included in the definition 
of the coherence function, the cross-spectrum bias error must 
be reduced to an acceptable level if a good statistical esti- 
mate of the coherence function is to be achieved. As stated 
above, the number of averages used in the impulse technique 
is usually not great enough to significantly reduce the bias 
error. However, the coherence function is still useful for 
indicating the importance of noise in the impulse technique. 
This is because noise in the signals causes variance in the 
value of the coherence function with frequency. This effect 
is illustrated in the section on measurement procedures. 

We oisplay Of Frequency Response 

The frequency response function is complex--that is, 
it has associated with it both magnitude and phase. There- 
fore, it can be displayed in a number of forms, including 
magnitude and phase versus frequency, real and imaginary 
Magnitudes versus frequency, and imaginary magnitude versus 
real magnitude. Each of these types of displays has its 
Own particular usefulness. The most common type of display 
for structural frequency response data iS magnitude and fre- 
quency plotted logarithmically. This type of display, with 
the magnitude in terms of compliance (ratio of displacement 
to force) is called a Bode plot. In this form of the fre- 
quency response function, resonances occur as peaks in com- 
pliance plots (points of maximum dynamic weakness) and all 


resonance peaks of equal damping have the same width regardless 
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of resonance frequency. Lines of constant dynamic stiffness 
have zero slope, and mass-dominated frequency response lines 
have a -12 dB-per-octave slope. Resonances occur as nearly 
circular arcs in the complex plane (real versus imaginary 
plot) with frequency increasing in a clockwise direction around 
the arc. In the case of real normal modes (which occur in 
systems with relatively low damping and with resonances well- 
separated in frequency), each resonance arc is approximately 
tangent with, and lies below, the real axis and is symmetric 
about the imaginary axis when the frequency response is ex- 
pressed as compliance. The complex plane plot is useful when 
certain types of analytical curve fitting operations are 

being performed on the frequency response data. The plots 

of the real and imaginary magnitudes of frequency response 
versus frequency are most useful when dealing with real normal 
modes. In this case the resonances will occur as peaks in 

the imaginary magnitude plot and the real magnitude will pass 
through zero at the resonance frequency when the frequency 
response is expressed as compliance. 

The frequency response characteristics of a structural 
element are determined by measuring a set of cross-frequency 
response functions as discussed in Reference 9. The cross- 
frequency response functions may be obtained by exciting at 
one location on the structure and measuring response at 
various locations, or by measuring the response at a single 


location to excitation at various locations. The resulting 
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frequency response functions comprise one column of the trans- 
fer matrix in the first case, and one row of the transfer 
matrix in the second case. Either set will, in general, com- 
pletely define the modal characteristics of the structural 
element. In mathematical terms the set of frequency response 
functions yields the eigenvalues and eigenvectors, which are, 
in general, complex terms. The real part of an eigenvalue is 
the damping and the imaginary part is the frequency associated 
with a given resonance. Each eigenvector defines a resonance 
mode shape. 

With real normal modes, each point on a structure is 
either exactly in-phase or exactly 180 degrees out-of-phase 
with any other point at the resonance frequency. Certain 
types of damping which are often encountered in practice will 
Cause the eigenvectors to have non-zero imaginary components, 
resulting in complex mode shapes. When a mode is complex, 
the relative phase associated with a point on a structure is 
some value other than 0 or 180 degrees, with the result that 
node lines (lines of zero deflection) are not stationary. 
Precise description of complex modes requires that some type 
mmanalytical curve fitting technique be applied to the fre- 
quency response data. 

The frequency response function of an operating system 
can be computed if the system input and output signals meet 
previously stated requirements of Fourier transformability 


and non-zero value, assuming the system input and response 
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can be measured. However, in practice there are usually 
multiple inputs to the system--either several inputs at 
different locations or inputs in more than one direction at 
a given location. In the case of multiple coherent inputs, 
the complexity of the analysis is greatly increased. For 
this reason, and the difficulty of accurately monitoring 
operating inputs, frequency response measurements are usually 
made by applying the system input "artificially" through 
some type of exciter. It is in the form of the input signal 
and the way it is applied to the structure that the wide 
variety of frequency response testing techniques arises. 

The usefulness of the impulse technique lies in the 
fact that the energy in an impulse is distributed continu- 
ously in the frequency domain rather than occurring at dis- 
crete spectral lines as in the case of periodic signals. 
Thus, an impulse force will excite all resonances within its 
useful frequency range. The extent of the useful frequency 
range of an impulse is a function of the shape of the im- 
pulse and its time duration. For a square pulse the frequen- 
cies of the zero crossings are at integral multiples of the 
inverse of the time duration of the impulse, illustrating 
the very important inverse relationship between the time 
duration of an impulse and its frequency content. 

The useful frequency range of an impulse is also a 
function of the shape of the impulse. By varying the weight 


and hardness of an impacting device and the manner in which 
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the impact 1S applied, the shape and time duration of the 
impulse produced can be varied to suit the measurement 
requirements. 

a. Nonlinearities in the Structures 

Excitation of a nonlinear system by a pure-random 
Signal will yield the best estimate (in a mean-square sense) 
of the linear system response. Excitation by a pure sine 
wave is also useful for studying nonlinear systems because 
it allows precise control of the input spectrum level. How- 
ever, the impulse technique, because of its very high ratio 
of peak level to total energy, is particularly ill-suited 
for testing nonlinear systems. Therefore, it is important 
to understand the various types of nonlinearities that can 
occur in structural systems and to be able to recognize 
nonlinearities in measured frequency response functions. 

One of the most common types of nonlinearities 
encountered in structures is that due to clearance between 
parts. This type of nonlinearity is frequently encountered, 
for example, when testing gear systems and shafts mounted in 
bearings. The effects of this type of nonlinearity on measured 
frequency response functions when using impulse excitation 
are poor estimates of static stiffness values and poor re- 
peatability of the frequency response estimates. Also, the 
apparent damping in the estimates will be greater than the 


actual examples. 
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The best method of dealing with this type of 
nonlinearity is to preload the system to take up clearances. 
Care must be taken when this is done, however, Seeeu ae any 
preload will change the boundary conditions of the structure 
and can itself lead to erroneous frequency response esti- 
mates. The usual approach is to apply the preload through a 
very soft spring so that the resonances associated with the 
meeltoad lie below the frequency range of interest. 

Another type of nonlinearity that is frequently 
encountered is nonlinear damping. Nonlinear damping effects 
are usually associated with joints in the structure, where 
the damping is a function of the relative displacement at the 
joint. In general, the frequency response estimates obtained 
by the impulse technique will agree most closely with those 
obtained with a low level of continuous excitation. However, 
if the point of excitation is close to a location where non- 
linear damping occurs, there will be high relative motion at 
that location, and the apparent damping in the measured fre- 
quency response will be high. In systems with low damping, 
this will give the measured frequency response a discontinuous 
appearance, due to the varying level of damping as the 
response to the impulse attenuates with time. 

The third type of nonlinearity that commonly 
occurs in structures is load-sensitive stiffness, where the 
Spring rate of elastic elements either increases or decreases 


meee toad. The most direct way to identify this type of 





nonlinearity is to measure frequency response as a function 


of static preload and observe the change in resonance 


frequencies. 


b. Signal Processing 


The particular characteristics of an impulsive 


force signal and the resulting 


structural response signal 


Make the impulse technique especially susceptible to two 


problems: 


hNetsowana Eruncaetton errors. 


While these problems 


occur to some extent with other frequency response testing 


techniques, their unique importance in the impulse technique 


requires special signal processing methods. 


itwveas pOlmLed Out 
the usable frequency range for 
shape and time duration of the 


that there is sufficient force 


in the previous section that 
an impulse depends on the 
impulse. In order to insure 


over the frequency range of 


interest, it is necessary that the first zero crossing of 
the Fourier transform of the impulse be well above the maxi- 
mum frequency of interest. For a given time duration the 

first zero crossing occurs at the lowest frequency for a square 
pulse. For that type of pulse the first zero crossing occurs 
at a frequency equal to the inverse of the time duration. 

A good rule of thumb, then, is to insure that the duration 
of the impulse is less than 2At, where t is the sampling 
interval in the analog-to-digital conversion process. This 
would put the first zero crossing of the Fourier transform 


of a square pulse at the Nyquist folding frequency, and the 
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first zero crossing of other pulse shapes above the Nyquist 
folding frequency. 

The sample length is equal to NAt where N is the 
number of digital values in each sample. A typical value of 
N is 1024. Thus, the duration of the impulse is very short 
relative to the sample length. This means that the total 
energy of noise represented in the time-sample can be on the 
order of the energy of the impulse, even for high signal-to- 
noise ratios. The noise problem is further aggravated when 
employing the zoom transform, which yields increased resolu- 
tion in a given frequency band by effectively increasing the 
sample length. 

With other techniques, the effects of noise are 
reduced by averaging the power spectrum and cross-spectrum 
functions prior to the computation of the frequency response 
function. However, only a few averages are usually used in 
miesimpulse technique. Otherwise, the time advantage of the 
technique is lost. Therefore, special time-sample windows 
have been developed for the impulse technique. 

At first thought it might seem appropriate to 
just set all time-sample values beyond the impulse to zero, 
Since it is known that the true signal value after the 
impulse is zero. However, this would be equivalent to multi- 
plying the signal by a narrow rectangular window. In applying 
any type of window, it 1S important to keep in mind that 


multiplication by a window in one domain is equivalent to 
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convolution of the Fourier transforms of the window and the 
data in the other domain, resulting in distortion of the 
transformed signal. This distortion will be mamumisea by 
minimizing the width of the main lobe of the window transform 
and suppressing its side lobes. However, there is a funda- 
mental conflict between these requirements and the reduction 
of noise in the time-sample because both the width of the 
main lobe and the amount of noise reduction are inversely 
proportional to the width of the window in the time domain. 
To further complicate the situation, suppression of the side 
lobes is generally achieved at the expense of broadening the 
main lobe. 

A good compromise has been arrived at in practice 
in the form of a window with unity amplitude for the duration 
of the impulse and a cosine taper, with a duration of 1/16 
Seetne sample time, from unity to zero. 

Noise problems may also be encountered in the 
response signal, particularly when dealing with heavily damped 
systems and when using zoom transform analysis. In both 
cases the duration of the response signal may be short rela- 
tive to the total sample time, so that noise may comprise a 
Significant portion of the total energy in the time-sample 
even with relatively high signal-to-noise ratios. Another 
error in the response signal that is encountered when testing 
lightly damped structures occurs when the response signal does 


not significantly decay in the sample window. In this case 
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the resulting time-sample is equivalent to multiplying the 
true response signal by a rectangular window, with the result 
that the frequency resolution may not be sufficient to resolve 
individual resonances. 

An exponential window has been developed to re- 
duce the errors that occur in both situations described above. 
The window decays exponentially from 1 to a value of 0.05 in 
the sample time. It can be applied directly to the time- 
sample of the response signal or to the impulse response 
function. As with all windows, the exponential window does 
change the resulting measured frequency response function; 
but its only effect is to increase the apparent damping in 
the resonances. It does not change the resonance frequencies 
and, because the effect of the exponential window is the 
same on all frequency response measurements, it will not alter 
the measured mode shapes if applied to all measured frequency 
response functions. In addition to reducing noise and trun- 
cation errors, the exponential window will also reduce errors 
which often occur when testing lightly damped systems in 
which the damping varies with the measurement position on 
the structure. 

Because the exponential window increases the 
apparent deamping in the resonance modes, there is a tendency 
of the window to couple closely spaced resonance modes. Zoom 
transform analysis may be required in some cases to allow 
sufficient resolution of closely spaced modes when using the 


exponential window. 
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Zoom transform analysis is discussed in some de- 
tail along with several examples in Reference 9. It isa 
very valuable tool in impulse testing, as it is in other 
frequency response measurement techniques. The effect of the 
zoom transform is to increase the resolution of the analysis 
by allowing independent selection of the upper and lower 
frequency limits of the analysis band. With the zoom trans- 
form, for example, it 1S possible to perform an analysis in 
the frequency range from 900 to 1000 Hz as opposed to the 
corresponding base-band range of 0 to 1000 Hz, resulting in 
a 10-to-l increase in resolution, for a given sample size N, 
in the 900 to 1000 Hz band. Because of greatly increased 
resolution possible with the zoom transform, it can be effec- 
tively used in frequency response testing to separate closely 
Spaced resonance modes. 

There are two important effects of the zoom 
transform in the impulse technique, both associated with the 
resulting increase in sample time. The first effect is to 
make possible much better estimates of damping in lightly 
damped systems. This is due to the reduction of truncation 
errors in the sampled response signal. The second effect, 
mentioned previously, iS aggravation of the noise problem 
in both the input and response signals. The second effect 
makes it essential that force and response windows be applied 
to the data in most cases when using the zoom transform with 


the impulse technique. 
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3. HP=-5451C Fourier Transfer Function 
A transfer function, as determined by the HP-5451C 
Fourier Analyzer, is a mathematical description of a system. 
It can be defined as: 


Fourier transform of output 


transfer function = 
hoUmrom thane tOorm OL Tnout 


or equivalently, 


average cross power spectrum 


Branster function = of input and output 


average power spectrum 
of input 
The coherence function measures the degree of causality 
between any two signals. It can, therefore, be used to check 
the validity of the transfer function. When a transfer func- 
tion is computed it may not be obvious that there are ex- 
traneous inputs, or that the system is nonlinear. Both of 
these factors would introduce error in the computed transfer 
function. The “Transfer Function" program of the HP-5451C 
1s used to compute the transfer and coherence functions. A 
program flow chart and a program listing are provided in 


Bepenaix B. 
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IV. THEORY OF THE TEST PROCEDURE TO MEASURE DAMPING 


The design concept of the experimental procedure is based 
on the Hewlett-Packard modal analysis software specifically 
designed for the HP-5451C Fourier Analyzer. The theory 
of the complex modes for damped oscillatory mechanical sys- 
tems is described below, much of which is explained in 
greater detail in the Modal Analysis Operating and Service 


Manual (Option 402). 


A. MODAL THEORY OF OPERATION [REF. 6] 
Assume that the motion of a physical system can be des- 
cribed by a set of n simultaneous second-order linear dif- 


ferential equations in the time domain, given by, 


Mx + Cx + Kx = f£ (9) 


where the dots denote differentiation with respect to time. 


f= yt 


1s the applied force vector, and 


an = Se) 


is the resulting displacement vector, while M, C, and K are 


the (nxn) mass, damping, and stiffness matrices respectively. 


oye, 





Our attention will be limited to symmetric matrices, and to 
real element values inM, C, and K. 


Taking the Laplace transform of the system equations gives 
Biegects) =) bis), (10) 


where: 
B(s) = Ms” + Cs + K (11) 


Here s is the Laplace variable, and F(s) is the applied force 
vector and X(s) is the resulting displacement vector in the 
Laplace domain. B(s) is called the system matrix, and the 


transfer matrix H({s) is defined as 


i 


Hc) — Es) e (12) 


which implies that 
Cs, his) = X(s) (13) 


Fach element of the transfer matrix is a transfer function. 
The elements of B are quadratic functions of s, and since 
H = Bt it follows that the elements of H are rational 
fractions ins, with det(B) as the denominator. Thus, H(s) 
can always be represented in partial fraction form. 

If it is assumed that the poles of H, i.e., the roots of 


det(B) = 0, are of unit multiplicity, then H can be expressed 


as 


ah 





’ 1) (14) 





The poles occur at s = P, (zeros of det(B)), and each pole 


k 
has n (nxn) residue matrix a, associated With at. For an 


aa order oscillatory system, there will always be 2n poles, 


but they will appear in complex conjugate pairs. Each com- 
plex pair of poles causes a mode of vibration in the struc- 


ture. The poles are complex numbers expressed as 


e =e + 1) E>) 


where OL is the damping coefficient (a negative number for 


stable systems), and w, is the natural frequency of oscilla- 


k 


tion. The resonant frequency is given by 


ae 2 
Qy = OL + wy (rad/sec) (16) 
and the damping factor is 
O 
C= ta ap 
k 
These coordinates are shown in Figure 7. When ¢ = 1, mode (k) 


is said to be critically damped. It is also possible that 
Cc > 1. For this case the poles of mode (k) lie along the 
real (or damping) axis in the S-plane and it is said to be 


Super-critically damped. 


a2 





(IMAG.) 






FREQUENCY 
AXIS 


Wy  (S-PLANE) 


(REAL) 


x 
CONJUGATE 
POLE 


Figure 7. Poles of a Mode (k) 


Modal vectors (uy) are now defined as solutions to the 


homogeneous equation 


where: 


By = B(p)) 


Bemeatll k = 1l,...,2n. 
In addition, pre-multiplying B times the equation (14) 


BeeeH(s), multiplying by the scalar (S-p,) and letting s = Py 


gives 
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(19) 


ft £ollows by comparison of (10) with (19) that every 


eolumn of ay must contain the modal vector u for each 


Kk! 
f= 1,...,2n-e 

Similarly, post multiplying B times the equation (14) 
fon H(s), multiplying by the scalar (S-py), and letting 


oe P), gives 


a, By = 0 (20) 
mereeeach k = 1,...,2n. This can be rewritten as 
oe 
By ay = 0 (249 


where t denotes the transpose. 


But Bue and hence Are is assumed to be symmetric so 
equation (21) is the same as (19) and it follows that each 
row of ay must also contain the vector u, HOt, Cac, kw = 1 eee ZN. 


In order to satisfy the above conditions a, must have the 


form 
a = Auut (n xn) (22) 
k Kak ek 
where A, is a scalar. 
This pervasiveness of the modal vectors throughout the 
transfer function matrix is evidence of the so-called global 
property of a mode of vibration. 
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In these terms, H can be rewritten as 


2n A 
H = ) — wu, (23) 
k=] 5° Py 





and this is easily written in matrix form as 


2S Come (mex 1) (24) 


where the columns of U comprise the u, modal vectors: 


k 


(n x 2n) (25) 


and L is a diagonal matrix containing all s dependence 





calor 


"= = | fon x 2n) . (26) 


=~ Pon 
Pre-multiplying H by One equation (13), can be written as 
(USUL) (UF) = (U®x) (27) 


so that ie transforms the spatial vectors F and X to vectors 
UF and u'x in modal coordinates. Similarly u‘uL is the modal 
representation of H. Since B(P, )u, = 0, it follows that 


* a3 
B(P,) uy = 0, so the modal vector associated with the conjugate 


pole (Py) is ay (the conjugate of u,)- 
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Thus, the above U matrix always contains conjugate pairs 
of modal vectors, and the L matrix always contains elements 
corresponding to conjugate pole pairs along its diagonal. 


ithg Us is defined as that (nxn) part of U associated with 


* 


DOsitive poles, then Uy will correspond to the negative poles. 


Similarly, L can be broken into two parts, L, comprising the 


il 
positive poles, and L. comprising the negative poles. It 
can then be represented 


= * x et 
H = U,L,U, + U;LQU, (28) 


or in partitioned form as 





a so 
= ‘ y ; 
s LU, ; Uy! au , Oj; (uy: UL, (29) 
ssocss 
0 ! L, 
Each of these sub-matrices is (nxn) and only L, and L, 
ree Lunctions of s. 
‘ae an 
ae ST Py 
L, = = ; Lo = 
A A 
n n 
= eon* 
SiR a 
(30) 


Then H can be written as: 
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te Ay k ORE a 


2 : 
SP, 





n 
H = : uu 3a) 


Ke. Kk S"P, 
Fach element of the H matrix has a different zero in the 
S-plane, depending upon the values of A, and U, at each 
Beant, but the poles of each element of H are common, and 
occur at s = Py and s = eye 
1. Identification of Modal Parameters 
Because of the form of the ay matrix, only one row or 
column of the transfer matrix need be measured and analyzed, 
Since all processes of meaSuring the transfer matrix, the 
unknown parameters in equation (14) (1.e., the complex values 
of Py and the complex values of the elements of one row or 
column or the residue matrix ay) are identified. 
Once one row or column of ay has been identified, it 
is then possible to construct the rest of the rows and columns 


For example, if the ei co Lumn sor ay is given by aig! 


~~ = (32) 


kq AUG kUK 
where Udy is the aco component of the modal vector. Since 
oak is a scalar, it is clear that the vector of modal 


residues a 1s proportional to the modal vector u for each 


kq 
wer, ... , on. 


ed 


Since A, is a scaling constant it can be assumed that 


either A) = l1 or the square root of utu = 1 without loss of 


Generality in the following derivation. 


=) 





Suppose A, = 1, then 


a = uy (33) 


and for the ge column or row of a 


k 
a eee Ol (34) 


h 


Hence the ae element of eg is 


2 
a = u 55 
kqa ( qk? (35) 
Since ek is a scalar, equation (34) can be rewritten 
a 
— AS (36) 
qk 
and substituting this back into (33) gives 
c 
a, ia 
a _kq"kq_ (27) 
k Ga 2 
qk 
or using equation (35) 
ic 
a, a 
a, = _kq_ kg (38) 
a 
kqq 


Hence the entire matrix ay and therefore, the entire 
transfer function matrix H(s) can be constructed once one row 
or column of residues ay has been identified as well as the 


pole locations pk, for each k = 1,...,2n. The residue a, _., 
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memeaeeene Griving point of the structure, i.e., the point 
where the structure is excited. 
7 impulse Response of Complex Modes 

It has been shown that a mode of vibration is repre- 
sented by a complex conjugate pair of poles and a complex 
conjugate pair of modal vectors in the transfer matrix. 
Hence for a single mode of vibration (k) the transfer matrix 
is written 
A 


,(S) = SP, AE S=p, (n xn) ey) 





It is convenient to remove a factor of 21 from the 


residue matrix, that is to define another residue matrix ry 


such that 


= 2.a (40) 


Then the transfer matrix 1S written as 


*k *k 
PS) ne eee ee (41) 
k 2i(s P,) 2i(s P,) 
Each component of this matrix exhibits the rectangular 
(or coincident--quadrature) from shown in Figure 8 for each 


valued rh: 


ae 





REAL PART 


H, (s) 








FREQUENCY 


IMAGINARY PART 


H, (s) 


fy 
FREQUENCY 


Figure 8. Transfer Function for a Single Mode of 
Vibration 


Taking the inverse Laplace transform of (41) gives 


1g joes ag Sas 
= k k~ _ -k k 
rel Ga) = Da. = En e GA7,) 
=) ne 
x(t) = izle sine(w,t + o,) (43) 
where a, = the angle of the complex residue r 


k k° 


Equation (43) is the impulse response of a complex 


mode of vibration. The mode shape is defined by the magnitude 


ae and the phase angle (oy) at Gach point on the structure. 
' ’ 
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The mode shape matrix is printed in the Modal System in 


terms of these magnitude and phase angles. 


3. Modal Mass, Damping and Stiffness and Scaled Mode 
Shapes 


Recall that the assumed symmetry of H(s) along with 
the global nature of mode shapes implies that the residue 
matrix for each mode is proportional to see If Qk = 2iA,, 


then equation (22) becomes 


mr. = Q,u,u (n xn) (44) 


Once a row or column of each residue matrix ry has 
been identified, the mode shape vectors UL can be calculated 
to within a constant of proportionality. This vector can 
be scaled large or small, if desired, by suitable choice of 
Qy- @riteria for choosing Q, will emerge from a discussion 


of modal mass, damping, and stiffness. 


Consider a single degree of freedom system represented 


by 
mx + cx + kx =f (45) 
where m, c, and k are scalars. The Laplace Transform gives 
ae + cs + k)X(s) = F(s) (46) 
so that 
i 
X(s) m 
H(s) = = (47) 
PRY 24 Se 4 
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If the damping factor is less than one, the polynomial 


in the denominator can be factored. 


i 
H(s) = ———“—~ (48) 
(So) UsSje 
where 
- a 2 k 
ae” ee a2) 
~ ; 2 
p= ti EH (50) 
Da eo + 1, Gi 0, tie 0 goplle 
and 
x 
p = -0 - lw (a2) 


To express Equation (48) in the partial fraction form 


Seeeaquation (39), solve for a and a”. 


L re 
CC) a ae 
(s-p) (s-p”) (s-p) aoe 
(5) 
iL : 
m _ a (s-p) 
= (ao 4-220. (54) 
(s-p) (S550 aga 
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DE en 
iL 
Ms * 
-21W aa 
i i 
7 mu) mw) 
H(s) = Silssor ” TN 


The residue r can be found from Equations 


ec. (57). 


Oe 


(3S), 


(56) 


i) 


(aie, 


(58) 


Now note that Equations (50), (51), and (58) give m, 


Speand k in terms of r, w, and o. 


1 
m = | 
rw 
Co = 925m 
a Wwe ~ 5G 


CS) 


(60) 


foi) 


For multiple degree of freedom systems, the analogous 


definitions are made for each mode k. 





m = i 
k Q,W) 
Cy = 20, mM, 


(62) 


(63) 





(64) 


Equations (44), (62), (63), and (64) are used to 
calculate modal mass, damping, and stiffness and scaled mode 
shapes in the Print Step. Qy Vopaketerary SO lt can be 
chosen to give unit mass, unit uu, etc. 

The modal mass, stiffness, and damping coefficients 
can be interpreted as being the masses, dampers, and springs 
of decoupled, single degree of freedom systems. These systems 
are equivalent to the original system under a change of 
coordinates. 

feeeericasurement Implications of Modal Theory 

A fundamental assumption of modal testing is that a 
mode of vibration can be excited from anywhere on an elastic 
structure, except of course along its node lines (zero points) 
Pier e it can't be excited at all. This is another way of 
Stating the result derived earlier (1i.e., that the same modal 
vector--scaled by a different component of itself--1is con- 
tained in every row and column of the transfer matrix). In 
addition, modal frequency and damping are constants which can 
be identified in any element of the transfer matrix, l.e., 
any transfer function taken from the structure. 

It is important to recognize that this global mode 
shape concept exists within certain spatial boundaries, beyond 
which vibrations will not readily propagate. If two linear 
systems are completely isolated, then a single composite mode 


including both systems is not meaningful. 
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Conversely, it is important to include enough spatial 
points in the measured data set to describe all of the vibra- 
tion modes of interest. If some region of a bounded system 
is not monitoried or excited, or if points are not chosen 
sufficiently close together, then some modes cannot be ade- 
quately represented. 

When a system is represented by the partial fraction 
form of its transfer matrix, a closed form solution for the 
displacement at any point, for any combination of modes, is 
readily obtainable using simple matrix-vector multiplication. 
This is particularly helpful when the response is only to a 
few modes of interest. 

The primary purpose for using modal coordinates to 
describe the dynamics of a linear vibration structure is that 
1t drastically reduces the amount of time and effort neces- 
sary to measure the dynamics in a laboratory. The essence of 
the modal concept is that once one row or column of the 
transfer function matrix has been determined, the entire matrix 
and hence the entire dynamics of the structure can be specified. 

In order to obtain valid modal results with the Modal 
Analysis system, several important assumptions of the modal 
theory described here must be satisfied. 

1. The structure must exhibit the behavior of a linear 
system. If you are not able to successfully curve 
fit equation (14) to measured transfer function 


data, then the behavior may not be linear. 
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The structure must exhibit the reciprocity of symmetry 
property. This can be verified by comparing the 
transfer function obtained from an excitation measure- 
ment at point A and response measurement at point B 
with the transfer function obtained from an excitation 
at point B and response at point A. This should be 
the same. 

The structure must have distinct pole pairs. Since 
the Modal Analysis system uses analytical expression 
foo tO pertemm the curve fitting, it cannot fit re- 
peated poles (repeated roots of det(B) = 0 because 
this expression was derived under the assumption that 
the poles are unique. This condition can be difficult 
to detect in one set of transfer function measurements. 
To insure that the poles are unique you should measure 
a different row or column of the transfer function 
matrix, and compare the mode shapes from the two sets 
of measurements. If they are not identical there is 

a strong possibility of repeated roots. 

There may be more than one modal vector corresponding 
to a given pole, due to certain kinds of boundary 
Sendititons On Ehe structure, 1.e., symmetric conditions 
along several axes. To insure that this is not 
happening, you should compare modal vectors obtained 
from two or more rows or columns of the transfer 


Mee ci x . 
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The transfer function H(s) is defined in terms of an 
input force and a response displacement. Response 
velocity or acceleration can be converted to displace- 
ment. But non-force Modal Analysis system will calcu- 
late, print, and display modal data (i.e., poles and 
residues) for such measurements, however. The excep- 
tions are the four mode shape scaling procedures that 
produce modal mass damping and stiffness. The system 
logica will reject an attempt to calculate these values, 


unless motion in response to force was measured. 
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V. PROCEDURE OF TEST CHAMBER CHARACTERIZATION 


A. PROCEDURE OF DATA COLLECTION 

After construction of the test chamber, preservation and 
painting, installation of cooling/heating coils, and instal- 
lation of all insulation, the following is the procedure used 
for final characterization of the test chamber (with front 
panel removed). 

1. Draw a 2.54 inch x 2.54 inch grid on entire inside 
surface of the test chamber. This grid 1S used to record 
location of impact or excitation, as well as location of 
response transducer (Figure 9). 

2. Re-install exciter and support framework; re-install 
specimen fixture with specimen in place; and re-install 
test chamber top. Drill 1/4 inch hole through specimen for 
attachment of exciter driver rod. Install exciter driver 
rod with Force Transducer in place (Figure 10). NOTE: 
Acoustic exciter and force transducer will not be used in 
this work. For impulse hammer testing disconnect exciter 
from specimen. 

3. Using accelerometer mounting wax, attach pick-up trans- 
ducer to desired location. Connect signal conditioners and 
wiring as instructed in the operating manuals. System flow 
chart is presented in Figure ll. Impulse hammer, signal 


conditioners and typical attachment of pickup transducer are 
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IMPULSE RESPONSE 

















PECKUe 
TRANSDUCER 
SN 7531 






0-500 LB 
IMPULSE HAMMER 
SN 269 












SIGNAL 
CONDITIONER 


SN 3491 


SIGNAL 
CONDITIONER 


SN 3468 


HP 54440A 
LOW PASS FILTERS 


SELECT CODES 1&2 








HP 5451C 
FOURIER 
ANALYZER 





HP 2648A HP 7245B 
GRAPHICS PLOTTER/ 
TERMINAL PRINTER 

Figure 11. Impulse hammer technique flow chart 


iia 















presented in Figures 12 and 13. All calibration data 

of sensors and mounting and instructions are included 

in Appendix C. 

Turn on the HP-5451C Fourier Analyzer (HP-2648A Graphics 
Terminal will automatically be energized). Insert 
"Modal" disk into disk drive. Energize disk drive. 

When "Disk Ready" light comes on "Boot" the system as 
outlined in volume 1 of System Operating Manual (Figures 
14 and 15). Activate user assignable keys fl and f2 

on the graphics terminal to "Graphics Mode" and "Alpha 
Mode" respectively. Procedures for key assignments are 
found in volume 5 of the Computer Operating Manual. 

From this point on, the procedure for taking impulse 
data is interactive. Locally generated computer pro- 
grams, which are designed to allow the data generated 

by impulse hammer technique to be stored and later 
processed by the off-line Modal software, are presented 
in Appendix D. What follows is a typical step-by-step 
impulse data taking session assuming pick-up transducer 
at location L051 and hammer impulse location at L066 


The following abbreviations will be used: 


C = computer Keyboard entry 

D = graphics terminal display 
G = graphics terminal entry 
R= user response 


iz 
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Typical attachment of pick-up transducer 
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Step Number 


ZO. 


mel. 


ne? . 


3. 


14. 


Normal "Boot" of modal with fl & f2 programmed. 

D: SET UP DATA NUMBER? (1 to 5). | 

R;C: 2 [ENTER] (want to preserve Demonstration data 
in set up #l, set ups #2 through #5 are available). 

Dron mmr 2 i, FINISH(3) 2 

R;C: 2 [ENTER] (must enter particular transducer 
specifications and general information for 
a particular test). 

D: LABEL? 

R;G: TANK CHARACTERIZATION [RETURN] (title of 
particular test). 

Dee Lene Se. 6 (FERST,LAST) . 

ne oho TURN] (you May modify one or any 
sequential series of lines of the set up data). 

D: Ll. NO. OF MEASUREMENTS? 

R;G: 16 [RETURN] (A measurement is one record of 
storage on the disk. One record is 1024 words 
in length, and a test 1S n measurements stored 
in n contiguous records on the disk). 

D: STARTING LOCATION OF MEAS. STORAGE (1 to 485). 

Roo Deke LURN| (Storage location 1 through 255 are 
on the removable disk, and storage locations 
266 to 485 are on the fixed disk). 

D: DATA BLOCK SIZE? (64 to 2048). 

R;G: 1024 [RETURN] (block size of 1024 will give 


maximum zoom capability). 


aa 





Mone Ds 4. POLAR TRANS. FUNCT. DISPLAY? LINEAR(1) or 
LOG (2) MAGNITUDE. 
iy. R-G: 2 [RETURN]. 
Moe Ds 5S. TRANSDUCER UNITS? ENGLISH(1), METRIC(2). 
fm. RG ) (RERURN |: 
20. D: TRANSDUCER SCALE FACTORS (VALUE,UNITS) 
1 = MV/G, 2 = VM/IPS, 3 = MV/IN, 4 = MV/LB 
A INEUt? 
21. R;G: 10.6,4 [RETURN] (values are obtained from 
calibration data of sensors). 
maze OD: 8B. RESPONSE? 
Boe R-Gs 10.12,4 (RETURN). 
24. Ds: G. AMPLIFIER GAINS? 
A. cE Wie 
25. R;G: 1 [RETURN] (enter transducer amplifier switch 
setting). 
mee DD: B. RESPONSE? 
meee R;G: 1 [RETURN] 
fore Ds LIST(1), EDIT(2), FINISH(3)? 
NOTE: At this point (l) may be entered and all data may 
be checked. If any corrections are required, enter (2) 


and use procedures outlined in steps 5 through 27. 


Zo). 


RG: 3° |RETURN). 


nor : 


The "set up" portion of the test is now complete. 


The next step is the "measure" portion. 


a0). 


IBYe 


SEE NUMBER? “(i €o: 5). 
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Bi. 


oz. 


eS 


34. 


a. 


BO. 
BM). 


BS . 


R;C: 2 [ENTER] (the computer is now in the "measure" 
mode) . 

Dz: TRANS(1), RANDOM(2), DISK(3), USER(4), 

FOURIER/ZOOM/GRAPHICS (5) ? 

R;C: 5 [ENTER] (the HP 5451-C is now ready to set 
up for an impulse test). 

R;C: SELECT BLOCK SIZE (1024 MAX FOR ZOOM) 

SET ADC FREQUENCY RANGE AS DESIRED 

SET TREGGER TO INTERNAL 

R;C: [JUMP] O [SPACE] 1 [ENTER] (this command calls 

maemusSCigmeis@cmsams !,50,51,52, and 59 for use). 

TRANSFER FUNCTION MEASUREMENTS 

ARE HP FILTERS INSTALLED? 

Pe (O = NO l = YES) 

R;C: 1 [ENTER] 


D: CUT-OFF FREQUENCY (CHANNEL NO)? 


NOTE: To eliminate unwanted or meaningless data beyond 


range of impulse hammer, enter channel number above which 


the computer will not sampel data. 


Bo . 
40. 
41. 


uz. 


EXAMPLE: 0-500 lb. hammer useful range is from 
0-6000 Hz, with block size 1024 selected 
channels above 307 are meaningless. 

Rees 307 [ENTER] 

D: IMPACT (1) OR RANDOM (2) EXCITATION? 

RC (ENTER 


D: NO. OF AVERAGES? 


Ve) 





moa k;C: 4 [ENTER] 
44. D: BASEBAND OR ZOOM MEASUREMENT? 
(OQ = BASEBAND 1 = ZOOM) 
45. R;C: O [ENTER] (first measurement must be baseband). 
46. D: BASEBAND IMPACT MEASUREMENT 
PRESS CONTINUE FOR MEASUREMENT 
NOTE: At this point trigger and analog-Digital Converter 
levels need to be set to do this. The following keg 
strodes should be done. 
47. R;C: SET REPEAT/SINGLE SWITCH TO REPEAT 
PRESS ANALOG IN 
SET ADC RANGE SWITCHES AND TRIGGER SLOPE AND LEVEL 
SET ADC TO SINGLE WHEN READY AND PRESS CONTINUE 
moe RCs [CONTINUE]. 
NOTE: Impact 4 times. The computer will beep after each 
data entry. After the computer processes the data the 
ready light will come on and the data in Block 0 will 
be displayed on the small screen of the computer. 
Processed data are stored in the following locations and 
can be displayed by the key strokes: 
[DISPLAY]0 [ENTER] Log transfer function 
[DISPLAY]1 [ENTER] Coherence 
([DISPLAY]2 [ENTER] Input power spectrum 
[DISPLAY]3 [ENTER] Output power spectrum 
[DISPLAY]4 [ENTER] Cross power spectrum 


(DISPLAY]5 [ENTER] Raw data 
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To obtain hard copy of any of the above outputs execute the 
following key commands [USER[ 5821 [SPACE] 35 [ENTER]. Then 
Srcer [USER] [PLOT) [ENTER]. To return plotting to graphics 
execute [USER] 5821 [SPACE] 6 [ENTER], and when the Graphics 
Terminal is in the Graphics mode (depress fl) and [USER] 
[PLOT] [ENTER] is executed all plotting to be sent to the 
graphics terminal. To continue, f2 must be depressed. After 
all preliminary graphing is completed on baseband data press 
[continue]. Typical plots of baseband data are presented in 
Figures 16 through 20. 
49. D: PRESS CONTINUE WHEN READY 
eee R;C: [CONTINUE] 
sl. D: SAVE DATA? 
(O = NO l = YES) 
me. R:C: O[(ENTER] 
53. D: MAKE ANOTHER MEASUREMENT? 
() = NO 1 = YES). 
mee 6 6OR;-C; §6«66L [ENTER] 
595. D: NO. OF AVERAGES? 
eo. R:C: R[ENTER] 
SB’. D: BASEBAND OR ZOOM MEASUREMENT? 
() = BASEBAND 1 = ZOOM) 
So. R-C: L{ENTER] 
29. D: 200M IMPACT MEASUREMENT 
MOVE CURSOR TO START FREQUENCY 


Povoom vALUL s(oWlLlTCH REGISTER 11) 
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NOTE: The desired start and end frequency for the zoom 


measurement is not entered with use of cursor. 


ao . 


ill . 


OZ. 


D: "TYPICAL" ZOOM IMPACT MEASUREMENT 
MOVE CURSOR TO START FREQUENCY 


PRESS | VALUE” (SWITCH REGISTER 11) 


CHANNEL = 49.0000 
PREQUBIGY, =—957.0999 HZ 
AMPLITUDE = ~-700.0026 E =-9.0000 


MOVE CURSOR, LO ENDSrREQUENCY 
PRESS "VALUE" 


134.0000 


D: CHANNEL 


Zo be. 0000 THz 


FREQUENCY 
AMPLITUDE = -381.7491 E =9, 0000 
ANALYZE OLD OR NEW DATA? 


i OLD (FROM THROUGHPUT FILE) 


2 NEW (WILL WRITE TO THROUGHPUT FILE) 


PRESS CONTINUE FOR MEASUREMENT 


Eee 2) (ENE 


NOTE: Whenever zoom is being done new data is required. 


eS 


R;C: [CONTINUE] 


NOTE: For zoom measurement the impact should be a series 


of rapid impacts with a varying interval between impacts. 


The trigger light will remain lit thorughout the data 


taking session. 


64. 


D: ANALYZING THROUGHPUT DATA "TYPICAL" 
CNER FREQ 1479 HZ/DIV 125.0 
DF: 000.2441406 BLOCKS LEFT 25 


ZOOM POWER 4 





NOTE: After data has been analyzed and desired plots 
made as described in step 43, ZOOM data should be stored 
for later use by the off-line modal software. A rec- 
commended data sheet is shown in Figure 21. It is 
important to log number of measurements and start loca- 
tion of measurement storage. 

ooe)6«6D:lhCU PRESS 6 CONTINUE WHEN READY 

NOTE: Generally it is desired to make multiple tests at 
different locations at the same frequency range so that a 
modal analysis can be done at a later time. 
fee RC: [CONTINUE] 
67. D: SAVE DATA? 

(0 = No 1 = YES) 
NOTE: The data taking session can be repeated by entering 
sO". 
Soe R:C: Il1[(ENTER] 
69. D: DATA STORAGE NO.? 
moe R;C: XXX [ENTER] 
NOTE: Data storage registers 0 through 500 are available 
for permanent storage of the data. 
71. OD: MAKE ANOTHER MEASUREMENT? 
(O = NO I = YES) 

NOTE: From this point the steps 57 through 71 will be 
repeated for as many zoom measurements as desired. Each 
zoom measurement will be stored in a next higher record 


number with the begin record number known. 
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— oo oe SE ee seem ee ee oe OE ee ae ee oe ee eee SN ee ose ae eee oe ee eee oe ee SP ee ee oe eee ee eee ee ee oe ee eee a ee 


DATA 


~ | I I i | ! I j J j ! | i I I 1 ! 
oo I ! I I | i I I J 1 | I 1 I ! | 
wa a | I | ! ! | | I J | | 1 j I I I J 
eo ey | 1 I ! | ! | I I | 1 1 I I I I ! 
oo }! 1 I | I ! j I J J 1 i ! ! I I i 
Ju | 1 I | I i I 1 I | | 1 | | ! I 1 
sso I I I [ J | i I | ! ! | I 1 I 1 

- | 1 I ! | J | I I ! J 1 | 1 | I ! 
te I I I 1 I i ! J I | 1 J I ! | ! 1 
igte <2 | I I I I I 1 | I ! | ' I | | ! 1 
oc - | i I | I | ! J I | J \ I i ! I | 
rear i 1 I ! I 1 I I ! j 1 1 I I J I I 
8) ~ | I ! ! I \ | I I | | I i I 1 I I 
eR ay ! I I | I 1 i I I ! i ' I ! ! I ! 
WZ ! 1 ! I I I ! J I I I I I i | i I 
- & ! | I ! i ] j I { | I 1 I i | I ! 
_ 0 1 1 I I I | ! I I 1 1 1 i I | I I 
et > I i 1 I I I I J I I ! I I I 1 I I 
i ! i I 1 I J \ I I I I I i i J I i 
~ | I [ I 1 1 j I I I | ' 1 I 1 1 i 

ul ae I i | 1 I j I I I i ! ! I | I I 
3) ao |! 1 i ! 1 I | I 1 ! j j I i | ! I 
z= ee 1 I I ! I \ ! I 1 I ! ' I I | I i 
ou J i I I I ! i I ! ! ! t I ! ! I ! 
wood ! i I I | I | | I ! 1 1 | I | 1 ! 
OaoH | i I I I I I I ! I I ' | | ! I \ 
ul i | I I I I | I I I I 1 I 1 i 1 I ! 
ce Oo | { I I I I I I I i I I I I I 1 j 
- | I l 1 1 I J I I | I ! 1 J I I ! 
on i I i 1 \ I I I I I I I i I ! ! 
Ce I I i I 1 I I I I I 1 I 1 i | I | 
ao ! I i i I I I I I I ! i | | | I ! 
ee ! i i ‘I I | ! j I I ! 1 I | | 1 ! 
| 1 i ! | I | ! I I ! ! I | | I I 

re Q | I I I I 1 I ] ! I ! \ I | ! I | 
aE | J ! I I 1 | ! I I i ! I | 1 I | 
ch. oc | i 1 ! I I I 1 I I i ! I I I I | 
re 8, I I I 1 ! I I 1 I ! I 1 I I ! 
Loa! bad I I I I I I I 1 I i I I I I | 1 ! 
~ I ! | i 1 1 I 1 ! 1 I I I | | I I 
GS I i I i I I I I 1 I I I I I I ! 
uJ Ld ~ I ! i I I I j I 1 I 1 I I I I ' 
a) tc Le I I i ! I I I ! ! ! ! I | I I ! 
zucoOa | I i I I i I I 1 I ! I I I I I 
xc -90 I I I I ! I | ! i i | I 1 I 1 I 
im >< u. | I I I I I J I I 1 { I I I I ! 
<x OO I I I I j I I 1 I I i I I I 1 l 
x al I I ! I I I ! I i I ! ! I I I I 
~ I I 1 I I I I I I I i J I I I 1 \ 
~ I I I 1 I I I I ! I I ! I I I I I 
lu N I I ! I I I I ! i I I I I I I I ! 
a t I I I I I I 1 I I I I I I I ! I I 
> od I 1 I ! | I ! j I I I j I I i I [ 
- I I ! I I | i | I I I | I i I i ! 
~ 1 1 ! 1 1 1 I ! 1 I I | I I I I I 
uJ 1 I I I I ! | I 1 i l J I i ! I I 
= I I ! i 1 | i | 1 I I J I I I I I 
< I I ! i i ! i | I | I ' I I I I I 
ao I I I I I I i i I i I J I I 1 I I 
=o I 1 I I I ! i ! i i 1 I 1 1 1 I I 
moO ( I I 1 I I I I 1 i I I I 1 1 1 1 
uj 1 1 i 1 I I | I 1 1 I | I 1 1 I I 
= | I I 1 I J I J 1 1 I i | 1 I 1 I 
o. I i i I | I 1 I i 1 I i I I I I I 
= es ! | i I ! I I | I | I 1 I I I I I 
re O | 1 I I I I I I 1 I I I J I l I 1 
i) <= I 1 1 I I ! I 1 I | I I I I I I | 
w” i 1 1 I 1 i I 1 1 1 i I I i 1 i I 


Sample data sheet 


Pele iiees 2 1 
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After all data has been taken and recorded, and the 
number of measurements and start record measurement has 
been correctly entered in the set-up table, the system is 
ready to proceed. 

2. R:C: [SELECT STEP] Modal Analysis Controller 
foe oD: STEP NUMBER? 
74. %™R;C: ([X] [ENTER] where X = set up 
data number il.e., 
1 for example 
eee: Lior (l), EDIT(2), FINISH(3)? 
Mo. R;C: 3{ENTER] assuming all data correct. 
From this point on, the procedure for the mode identifi- 
cation and damping measurement is covered in Reference 6. 
For the tank characterization the modes identified and damping 


measured at the particular modes are presented in Table II. 


B. TEST CHAMBER MAJOR MODES AND DAMPING 

To identify the major modes and measure the damping of 
the test chamber the polar plot of the baseband band was 
divided into 10 sections (Figure 22), and each section was 
zoomed using the same pickup location and using the same im- 


pact location on left, right, deck and back panel. 
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TABLE I 


TEST CHAMBER FREQUENCY SECTIONS 


SECTION START FREQ. STOP FREQ. CNTR. FREG2 DELTA FREQ. 
NO. eras) (Hz) Hz) (Hz) 
i 423 ue lowe! Teal 30) 5 
Zz 976 ea o llearea7 2a I Or ora 
5 ie oS 2,297 1,780 1.786 
4 PEN a 2S 2,481 Me a 7-al 
5 2a oD Sch B) Os oie 
6 3,540 Oil: 3520 0) a 
7 Ser bl 4,493 aly, 10) 73 O53 
8 4,473 4,805 4,639 0.814 
9 4,805 5,411 57 Lg ime630 
TO 5,411 3 vas Sy ayes 7 0.814 


A modal analysis of each section was completed. Major 
modes and damping for each major mode was analyzed. The 
results are presented in Table II. Test data are included 


in Appendix E. 
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SECTION 
NO. 


10 


NO Cor 
MODES 


ar2 


10 


10 


FABLE IT 
TEST CHAMBER MODES AND DAMPING FACTORS 


AVERAGE DAMPING 
FACTORS (3%) 


jaye sFe) 
RO ok 
woo re 
e545 
BS 
ge Tas. 
Pon 
4s ALS, 
leon 


meEGGo 


STANDARD DEVIATION 
OF DAMPING 


AsO) pis, 
OOK 
.2458 
.5084 
Oe 
5 0 ts 
BAA Te 
Woz 7 
Otog 


nv aZS 


The overall average damping factor of the test chamber is 


0.46673 with an overall standard deviation of 0.36893. 
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VI. PROCEDURE FOR DAMPING MEASUREMENT OF SPECIMEN 


pee SPECIMEN SECURED IN FIXTURE INSIDE TEST CHAMBER 

The test specimen (cast nickel-aluminum bronze, code 
FTC) was secured in the test chamber, impact locations A, B, 
C, are shown in Figure 23. To identify the major modes and 
measure the damping factor of the specimen the 0 to 500 lb. 
impulse hammer was used for the frequency range up to approxi- 
mately 6000 Hz, and the 0 to 50 lb. impulse hammer was used 
for the frequency range from approximately 5000 Hz to approxi- 
mately 12,000 Hz. Baseband polar plots for each impulse 
hammer at each location are shown in Figures 24 through 29. 

It should be noted that there is very good correlation between 
the 0-500 lb. impulse hammer and the 0-50 lb. impulse hammer. 
Additionally, there is very good correlation between various 
maeact locations. 

Utilizing the baseband data (0-500 lb. impulse hammer), 
obtained for impact location B, five frequency ranges were 
investigated. Additionally, using the baseband data obtained 
from the 0-50 lb. impulse hammer at location B, the frequency 
ranges of three additional ranges were investigated. Specimen 
Frequency sections are listed in Table III. Specimen modes 
and damping factors are presented in Table IV. Measurement 


data and rough log sheets are included in Appendix F. 
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TABLE III 


SPECIMEN FREQUENCY SECTIONS (SPECIMEN FIXED) 


SECTION START FREQ. SLOP GEREO. CNTR. PREO; DELLAS? REQ. 
NO. (Hz) (Hz) (Hz) (Hz) 
ik 410 Sho 673 IEP SSS as, 
2 840 Ee BS) 8 JE 7 Sie 4.882 
3 oA eye) s: 3,340 J 2) Sie) eS Bo) 
4 4,043 4,512 4,277 ieee 2 
5 4,747 27,0208 Se lee lh 9 53 
6 Dua (es 6,494 4.0469 
7 i271 o So) le 210 8,448 6.104 
8 Seco ioc OF 200 4.883 


Modal analysis of the specimen sections was completed and 


the results are shown in the following table 


HOTZ 





TABLE IV 


SPECIMEN MODES AND DAMPING FACTORS (SPECIMEN FIXED) 


SECTION NO. OF AVERAGE DAMPING STANDARD DEVIATION 
NO. MODES FACTORS (3%) OF DAMPING FACTOR 
It 1 O75 5 ae 
2 2 Oz /3 0.1088 
5 5 0.0559 Ols 0063 
= 2 0.0883 0. 00 25> 
5 2 Ce) Sie OR O576 
6 6 OL hehe. 0.0542 
7 5 050917 0.0541 
8 5 0.1000 Ors0S9'6 


The overall average damping factor for the specimen is 
about 0.1112% with an overall standard deviation of about 


0.06013. 
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B. SPECIMEN REMOVED FROM FIXTURE AND TEST CHAMBER 

In order to ascertain the effects of the specimen fix- 
ture and test chamber, one additional test was conducted on 
the same specimen. During this test the specimen was removed 
from the fixture and test chamber. The specimen was laid 
flat on 3/4 inch foam rubber. No other support was provided. 

Response and impulse hammer locations are identical to 
the previous test. Frequency ranges for use of the large 
(0-500 lb.) impulse hammer and small (0-50 1b.) impulse hammer 
remained unchanged. 

The baseband plots for large and small impulse hammer 
for locations A, B, and C are presented in Figures 30 through 
35. For the large impulse hammer four ranges of zoom were 
completed, and for the small impulse hammer five additional 
sections were zoomed. Specimen frequency sections are tested 


in Table V. 
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TABLE V 


SPECIMEN FREQUENCY SECTIONS (SPECIMEN FREE) 


SECTION START FREQ. STOP FREQ. CNTR FREQ. DELTA FREQ. 
NO. eis) (Hz) (Hz) (Hz) 
il 1094 1954 WA 2e! Ds WON al 
2 OES 8555 2755 4.883 
3 Sa ies 4610 4092 Seo pe 
4 4610 5293 4951 053 
5 4688 6299 5493 4.069 
6 6250 V2 2g, 6738 Dey) 
7 7276 8594 7935 Seaa 
8 8545 9473 9009 Daa 
9 9424 11,190 Ors 07 4.883 


Modal analysis was compled on the free specimen, and the 


results are tabulated below in Table VI. 
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TABLE VI 


SPECIMEN MODES AND DAMPING FACTORS (SPECIMEN FREE) 


SECTION NO. OF AVERAGE DAMPING STANDARD DEVIATION 

NO. MODES FACTORS (3) OF DAMPING FACTOR 
i 2 0.0885 0.0064 
2 5 .1264 nO 74 
3 1 .0339 N/A 
4 3 .0265 0017 
5 2 .0295 .0020 
6 4 .0250 .0080 
7 5 hes .0195 
8 3 0209 0123 
9 3 G0 se: 0.0069 


Again, all raw data sheets are included in Appendix F. 
The overall average damping factor for the free specimen is 
about 0.0500% with an overall standard deviation of about 


0.05043. 
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VII. RESULTS AND CONCLUSIONS 


Mode identification was accomplished utilizing the zoom 
feature of the HP-5451C Fourier Analyzer. After a baseband 
test was conducted and results stored (Figure 22), any section 
could be chosen to zoom. For example, zooming section 1 of 
the test chamber results ina plot as shown in Figure 36. 

The initial display will be in polar form and must be con- 
verted to rectangular (one button function on the HP-5451C) 

as shown in Figure 37. As shown in Figure 38, the coherence 
for this test was very good, and the data may be assumed to 

be accurate. For mode identification it is necessary to get 

a Nyquist plot of the data (in rectangular form). The Nyquist 
plot is vector representation of the real vs. imaginary values 
of the data. On a Nyquist a pure mode would sweep out a 
perfect circle. Modes with differing damping factors exhibit 
Circles with varying radius. Figure 39 shows the Nyquist 

plot of the zoom data. Note there are many modes present 

with greatly varying damping factors. Individual modes can 
Bemisolated by inspection of a small portion of the rectangu- 
lar data. For example, Figure 40 shows the rectangular plot 
meeme about 675 Hz to about 850 Hz. If the isolated portion 

is converted to Nyquist, as shown in Figure 41, the major mode 
at 772 HZ is readily identified. The Nyquist plots can be 
smoothed (made more circular) by higher resolution. That is, 


complete a zoom measurement on a smaller frequency band. 
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All modes identified for the test chamber and the specimen 
(both fixed and free) are listed in the tables included in 
Appendices E and F. 

It was noted that there are many more modes for the test 
chamber than for the specimen, but that is expected because 
the test chamber is a much more complex system, and the 
boundary conditions are considerably more complex. 

The number of modes for the specimen in the test chamber 
is found to be 26, whereas the number of modes of the speci- 
men separated from the test chamber is 28. This is good 
agreement, and it should be noted that the modes of the test 
chamber did not increase the number of observed modes of the 
specimen. 

The average damping factor observed for the test chamber 
memabout 0.47%, which is a typical value for plain carbon 
steel. The average damping factor for the specimen (fixed) 
is about 0.11%, and for the specimen (free) is about 0.05%. 
As expected, the fixing of one end of the specimen in the 
test chamber did increase the measured damping factor, but 
with more testing it is possible to quantify the increase 
and make suitable corrections. 

Perot Of Adamping factors vs. frequency for the test 
chamber is presented in Figure 42. The damping factor is 
higher at lower frequencies with the maximum damping factors 
being in the range of 1500-3000 Hz. After 3000 Hz the damping 


factor averages out to about 0.2% at higher frequencies. 
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For the specimen, damping factors vs. frequency (for 
both fixed and free) are presented in Figures 43 and 44. 
Like the test chamber, the damping factor appears to be 
higher at lower frequency for the fixed specimen, and 
without question this same trend is observed for the free 
specimen. The fixed specimen damping factor levels out to 
about 0.10% at higher frequencies, and for the free specimen 


it levels out to about 0.025%. 
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VIII. RECOMMENDATIONS FOR FUTURE WORK 


1. Expand the test procedure to include random excita- 
tion of test specimen for the energy input and measurement 
of responses up to 15,000 Hz. 

2. Conduct similar tests on various materials including 
high damping materials such as Incramute and Sonoston. 

3. Investigate the temperature dependence of viscous 
damping, especially for constrained layer and non-metallic 
composites. 

4. Investigate the changes in damping for materials 


immersed in a water medium at varying temperatures. 


12> 





APPENDIX A 


THE HP-5451C FOURIER ANALYZER 


The HP-5451C Fourier Analyzer System performs analysis 
of time and frequency data containing frequencies from DC 
to 50 KHz. The system analyzes time-series data such as 
mechanical vibrations, sonar echoes, tidal action, biomedi- 
cal phenomena such as brain waves and nerve impulses, 
voltages and currents in electronic systems, and acoustic 
phenomena. These analyses may detect signals hidden in noise, 
or may locate critical functions in complex systems. Both 
continuous and transient data may be processed. 

Keyboard programming allows the following operations 
automatically without special software: 

A. Forward and inverse Fourier transform 

Bs. Magnitude and phase spectrum 

C. Power and cross power spectrum 

DBD. Transfer function 

E. Coherence function 

me Convolution 

G. Auto and cross correlation 

H. Hanning and other weighting functions 

im Histogram 

wee Scaling 


K. Ensemble averaging (time and frequency) 
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Six editing keys operate an on-line resident editor so 
that a sequence of steps configured into an automatic 
measurement procedure may be changed on-line without ene 
need to do off-line editing, compiling and testing. [In 
fact, the series of steps or programs used to perform a 
particular operation can be stored on the Disc for easy 
re-entry into the Fourier Analyzer. 

Data input and output is likewise controlled from the 
keyboard. Data can be entered in analog form through the 
four channel Analog-to-Digital Converter, or in digital form. 

Results of all operations are displayed on the oscillo- 
scope. In addition, results can be printed out in decimal 
numbers or plotted on the Graphics Terminal or an X-Y plotter. 

The Fourier Analyzer is a completely calibrated system; 
all displays and data output are accompanied by a scale 
factor relating them to physical units. This calibration 


results from digital techniques being used in all computations. 
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HP-5451C SPECIFICATIONS AND CHARACTERISTICS 


SPECIFICATIONS 
(Specifications describe the standard system’s warranted performance.) 


ANALOG-TO-DIGITAL CONVERTER 


Input Voltage Range: +0.125V to +8V peak in steps of 2. 

Input Coupling: dc or ac. 

Input Channels: 2 channels wired for 4 standard. 4 channels op- 
tional with plug-in cards. 

Resolution: 12 bits including sign. 

input Frequency Range: dc to 50 kHz, 5 Hz to 50 kHz, ac 
coupled (100 kHz optional). 

Sample Rate: 
Internal: 100 kHz max. (1, 2, 3, or 4 channels simultaneously). 

(200 kHz optional on 1, 2, 3, of 4 channels.) 

(50 kHz max. (3 or 4 channels simultaneously).f] 

An external time base may be used to allow external 

control of the sampling rate up to 100 kHz (200 kHz 

optional). One sample can be taken for each clock 

pulse (TTL level). 

Internal Clock Accuracy: +0.01%. 


External: 


DISPLAY UNIT 


Vertical Scale Calibration: Data in memory is automatically scaled 
to give amaximum on-screen calibrated display. The scale factor 
is given in volts/division, volts2/division, or in dB offset. 

Log Display Range: 80 dB with a scale factor ranging from 0 to 
+998 dB. Offset selectable in 4 dB steps. 

Linear Display Range: +4 divisions with scale factor ranging 
from 1 X 10-512 to 5 X 10512 in steps of 1, 2, and 5. 

Digital UP/DOWN Scale: Allows 8 up-scale and 2 down-scale 
steps (calibrated continuous scale factor). 

Horizontal Scale Calibration: 

Linear Sweep Length: 10, 10.24 or 12.8 divisions. 
Log Horizontal: 0.5 decades/division. 
Markers: Intensity markers every 8th or every 32nd point. 


BASE SOFTWARE 


Transform Accuracy: The expected rms value of computational 
error introduced in either the forward or inverse FFT will not 
exceed 0.1% of the rms value of the transform result. 

Dynamic Range: >75 dB for a minimum detectable spectral com- 
ponent in the presence of one full scale spectral component 
after twenty ensemble averages for a block size of 1024. 


EXECUTION TIMES® 


Fourier Transform: <55 ms 
Stable Power Spectrum Average: <80 ms 
Stable Trt-Spectrum Average: <220 ms 


REAL TIME BANDWIDTHS* 


Fourier Transform: = >7.5 kHz 
Stable Power Spectrum Average: 5.4 kHz 
Stable Tri-Spectrum Average: 1.9 kHz 


MASS STORAGE SOFTWARE 
MAXIMUM REAL TIME DATA ACQUISITION RATE 
(Single Channel): 
BS 256: 10kHz 
8S 1024: 39 kHz 
8S 4096: 80 kHz 


(25 kHzft) 
(30 kHzf) 


OFF-LINE BSFA SOFTWARE 


Center Frequency Range: dc to one-half the Real Time Data 
Acquisition Rate. 

Center Frequency Resolution: Continuous resolution to the limit 
of the frequency accuracy for center frequencies >0.02% of the 
sampling frequency. 

Frequency Accuracy: 

Bandwidth Selection: {n steps of f/5n where n = 2, 3, 4, etc. 

Max. Resolution Enhancement: >400 

Dynamic Range:** 90 dB trom peak out-of-band spectral com- 
ponent to the peak level of the passband noise. 

80 dB from peak in-band spectral component to the peak level 
of the passband noise. 

Out-of-Band Rejection: >90 cB 

Passband Flatness of the Digital Filter: 


+0.01% 


+0.01 dB 


ENVIRONMENTAL CONDITIONS 
Temperature Range: 0°C to 40°C (104°F). 


*for band limited random noise type signals at block size 1024, no display, no Hanning. 
**After erght ensemble averages of a power spectrum at block size 1024 Reduced by 
10 dB at the exact center of the band. 
t These rates apply to systems with modules 54668 and 54451A/8 having a serial prefix 
lower than 1842. 
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SUPPLEMENTAL CHARACTERISTICS 
(Supplemental Characteristics are intended to provide useful information for 
system applications by giving typical, but not warranted, performance parameters.) 


ANALOG-TO-DIGITAL CONVERTER 


input Impedance: 1 M2M in parallel with <75 pf. 
Sample Rate Control: 

Maximum Frequency Mode: Maximum frequency selectable. 
from 0.1 Hz to 50 kHz (100 kHz optional) in steps of 1, 2.5, 5. 
This mode automatically sets maximum frequency indepen- 
dent of block size. 

Frequency Resolution Mode: Frequency resolution selectable 
from 0.2 mHz to 1000 Hz in steps of 1, 2, 5. This mode auto- 
matically sets frequency resolution and sample record length 
independent of block size. 

input Mode: There is a buffered and non-buffered analog mode. 

In the buffered mode, other operations can be performed on 

previously collected data while the ADC collects current input 

data into a buffer. 


DISPLAY UNIT 


Data may be displayed in single sweeps or refreshed continu- 
ously in the following forms: 


Amplitude 
Real Part 


Time (Linear or Log) 
Frequency (Linear or Log) 
Frequency (Linear or Log) 
Frequency ‘Linear or Log) 
Frequency ‘Linear or Log) 
Real Part 


- Imaginary Part 
Magnitude (Linear or Log) 
Phase 
Imaginary Part (NYQUIST PLOT) 





CRT Positioning: Three markers to aid in adjusting trace posi- 
tion as well as vertical and norizontal controls are provided for 
display positioning. 

Origin: Left edge of display, zero amplitude. 
+FS: Positive full scale, center of display. 
-FS: Negative full scale, center of display. 

Analog Plotter Output: Any displayed data can be output to a 
plotter or remote oscilloscope. HP 10640B interface 15 required 
for operation. 
Amplitude: 0.5V per oscilloscope display division. 
Linearity: 0.1% full scale. 
Interpolation: Linear interpolation in 0.05% steps. 

Type of Display: Points, bars, or continuous (interpolation). 


BASE SOFTWARE 


System Accuracy and Range: The Fourier Transform is imple- 
mented using conditional scaling for maximum accuracy with 
no data overflows allowed. All calculations use floating point 
arithmetic on a block basis with full 16- and 32-bit arithmetic 
where applicable. : 

Data Word Size: 16-bit imaginary with 32 bits preserved for 
double precision functions. Division, addition, or subtraction 
operations performed in 16 or 32 bits depending on data. 

Maximum Block Size: 4096 time domain points. 

Minimum Block Size: 64 time domain points. 

Permanent DataSpace: 28K words (16K words standard in systems 
with serial prefix below 1842, optionally expandable to 28K with 
option 011). 

Permanent Program Space: 


BSFA SOFTWARE 


Maximum BSFA Blocksize: 1024 time domain points. (2048 with 
option 670.) See also Table 5-1 in BSFA section of manual. 


32K words. 


MASS STORAGE 


Dise Unit: 


_ Capacity: = 2.45 megawords 
Data Transfer: 2.5 million bits/second 
Dises: 2 (1 fixed, 1 removable} 


% of Real Time at 100 kHz ADC Sampling Rate (Single Channel): 


BS 256: 10% BS 1024: 39% (25%T) BS 41096: 80% (30%) 
Number of Records Per File: 
Data Block: 214 records (4096 blocksize maximum/record). 
ADC Throughput: 199 records (4096 blocksize max./record). 
Program Stack: 138 records (470 steps/record). 
ASCII Text: 690 records (128 words/record). 
Index: 69 records (10 pointers/record). 
System Coreload: 4 records (32K words/record). 
Common: 286 records (256 words/record). 
Overlay: 20 overlays (8K words maximum, 7K words maximum 
with option 261 or 265). 


OPTIONS 261 & 265 
The magnetic tape options are used for ADC Throughput only. 


Maximum Real Time Data Acquisition Rate (Single Channel): 


Opt. 261 Opt. 265 
BS 256: 6kHz BS 256: 9kHz 
BS 1024: 12 kHz BS 1024: 21 kHz 
BS 4096: 15 kHz BS 4096: 30 kHz 


Number of Tracks: 9 

Read/Write Speed: 45 ips 

Density: Option 261, 800 bpi; Option 265, 1600 bpi. 

Data Transfer Rate: Option 261, 36K cps max.,; Option 265, 72K 
cps max. 

Rack Height: 610 mm (24 in.) 


POWER REQUIREMENTS, SIZE, WEIGHT 


Power Source: 115/230 volts +10%, 50/60 Hz. 1800 watts typical 
for base system. 
Size: Dimensions are for a typical system (excluding cabinet and 
terminal}. 
Height: 
Depth: 
Cabinet: 
Panel Height: 1422 mm (56 in.; 
Overall: 1631 mm (641% in.), height; 533 mm (27 in.}, width; 
762 mm (30 in.), depth. 
Weight: Net weight for a base system (excluding terminal) — 
163.3 kg (358 Ibs.). 
Price and Ordering Information: Consult the 5451C Fourier 
Analyzer System’s Ordering Intormation Guide. 


771 mm (28 in.); Width: 425 mm (1634 in.); 


616 mm (24V4 in.). 


SYSTEM INSTALLATION 
Included in the 5451C System is on-site installation. On installation, 
a trained Hewlett-Packard representative will perform an opera- 
tional demonstration to ensure that the system is functioning normally. 


TRAINING 
A course on Fourier analysis and system operation is optionally 
available at HP’s Santa Clara, California facility. On-site training 
can also be provided, if desired. 


tThese percentages apply to systems with modules 54668 and 544S1A/8 having a serial 
prefix lower than 1842. 
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APPENDIX B 


HP-5451C TRANSFER FUNCTION AND POWER SPECTRUM 
FLOW CHART AND PROGRAM LISTING 


This appendix contains flow charts and listings of the pre-written soft key programs shipped on the 5451C 
operating disc pack. These programs are for the Gold Key F2 (Transfer Function) and F5 (Power Spectrum) 
soft keys. The purposes for including these flow charts and listings are as follows: 


1. For those who wish to modify the programs to better fit specific applications. 
-2. As an example for those who wish to write their own soft key programs. 


3. As a model for re-entering portions that may have been accidentally written over on the disc 
(programs reside in unprotected areas). 


As an aid toward rapid understanding of these programs, a certain programming “style” has been used. As 
a result of this style, the programs are longer than they would otherwise need to be, but are easier to 
comprehend. Some elements of the style are described below. 


BLOCK STRUCTURE THROUGH USE OF LABELS 


The code in the programs is organized in “‘blocks” — functional segments which are delimited by LABEL 
instructions. A convention is followed in the choice of label numbers. The beginning of each block is 
designated by a label number ending in a multiple of 50, e.g. L1050. The block ends with the label 9 higher, 
e.g. L1059. Within the block, label numbers are between these limits, e.g. L105!, L1052 etc. One block may 
contain others — each will use the above delimiting convention. 


Although the use of these labels lengthen the programs, it makes them much easier to understand, and to 
correlate with the flow diagrams. For example, when several branches of a program go to the same point, 
each branch will jump to its own appropriate label. Since these labels all appear together at one point in 
the program, one can tell that multiple branches converge at that point. An example is in the Transfer 
Function program, stack 54, where labels 1209, 1259 and 1609 all appear together. This is just one example 
of the way in which complex internal linkages in the program are made more visible. 


BRANCHING THROUGH USE OF “COMPUTED GOTO’s” 


In most complex programs, branching is common. One means of branching is to use an IF statement, 
provided in Keyboard language by the “GOLD KEY” “SKIP” instruction. When there are more than two 
possible branches, however, use of IF branching tends to get complicated, involving multiple decision 
points. The “Computed GoTo” or “switch” type of branching statement is more suitable in such cases for 
simplicity of understanding. It has been used extensively in these keyboard programs, even for simple two 
way branching. By standardizing on it, the code becomes recognizable and easier to read. 
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In Keyboard language, the “Computed GoTo” is implemented by computing the number of a label and 
jumping to it. The following is an example. 


L 4050 Start of branch block 
YAt 0 14050 4D Set variable J to 1050 
+ (the value of variable 
44) 
J 0D Jump to the label number 
in variable 0 
L 4054 Code to be executed if 
j variable 441 = { 
J £0S9 Go to end of block 


beit52 Code to be executed if 
variable if =¢ 


L 1959 End of oranch block 


USE OF SUBROUTINES 


In these programs there are several functions which have been set up as subroutines. These inciude the 
parameter entry routine (see below) and the measurement routines. The: measurement routines are 
handled this way to simplify the flow of the main program and allow easy replacement of the 
measurement code in case you modify it in some way. 


PARAMETER ENTRY ROUTINE 


Since these programs ask the operator for many input parameters, a single subroutine, LABEL 100 in 
STACK 0, handles all the parameter entries. The routine is called with variable parameters 1 and 2 equal to 
the lower and upper limits on the range of allowable operator inputs. The routine reads your input, 
checks it against the range limits, and, if it is valid, passes it back to the calling program in (floating point) 
variable parameter 2000. If the input is out of range, this routine notifies you and waits for the new input. 
The routine will not return to the calling program until a valid input has been received. 


PRECAUTIONARY NOTES 
The following precautions apply to the operation of the preprogrammed measurements: 

1, When using the standard software zoom (BSFA), the measurement blocksize can be no larger than 
1024. When using the Option 670 Feurier Preprocessor for BSFA measurements, the maximum 
blocksize is 2048. 

2. The maximum center frequency you may specify for a BSFA measurement is 32767 Hz. 

3. The messages CF WHAT? or BW WHAT? may result if the center frequency and/or bandwidth you 


have chosen for your measurement are such that the BSFA analysis band is either less than 0 or 
greater than the ADC Fmax setting. Specifying different parameters shouid remove this problem. 
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4. The message DL WHAT? may occur when performing the on-line BSFA measurment. This is 
because the display is active during the on-line measurement. To remove this problem, either 
reduce the measurement bandwidth (thereby increasing the zoom power and lowering the data 
rate into the computer), or edit the appropriate keyboard stacks (stack 56 for transfer function, 
stack 61 for power spectrum) to remove parameter n3 from the calls to User Prog 45 far the on-line 
measurement (refer to commented program listings which follow). 


5. When performing an off-line BSFA measurement with an optional Mag Tape unit, perform the 
following steps before making the measurement. 


Set: 


ADC SAMPLE MODE to INTERNAL KH2z/ys 
MULTIPLIER to 100/10/5 

INPUT SELECTOR to A 

TRIGGERING to FREE RUN 

OVERLOAD VOLTAGE A to CHECK 


Enter: 


BLOCKSIZE 4096 ENTER 
MASS STORE 32 ENTER 
MASS STORE 22 SPACE 1 SPACE 150 ENTER 
MASS STORE 32 ENTER 


This writes 150 records of data on the magtape so that the magtape will be able to position to 
record 135 on the tape when the ADC throughput is performed. It will position by looking for the 
interrecord gaps written by the WRITE ADC throughput command. 


6. After completing a BSFA measurement, be sure that all data space declared by the zoom programs 
is released by pressing RESTART. 


As you go through the flow charts and commented listings, remember that these are only examples of 
programming the soft keys F1 through F6 on the Keyboard. It is up to you to determine which, if any, 
portions of these programs should be maintained. Because these programs are stored in unprotected 
areas of the Disc, there is the possibility they can be written over. If this should happen, you should enter 
the program stacks from the listing, substitute your own program, or copy from your back-up disc. 


The soft key programs and the associated ASCII text and variable parameters were originally stored on the 
system disc pack in Files 3,4, and 7. The records used are as follows: 


File 3 (Keyboard Programs) File 4 (Text Buffers) 
Record 0 Text buffers 51 through 55 
Records 51 through 62 . ASCII records 324Sthrough3444* 


File 7 (Common) 
Common) Record 0 
*This assumes that there are}70O records in File 4. If not, the first and last ASCII text records may be 
computed as follows: 
First record number = NR - (5 X last text buffer number) 
Last record number = NR - (5 X first text buffer number) + 4 


where NR = number of records in File 4. 


The allocations above should be kept in mind so.the above records are not inadvertently altered or 
destroyed when using the Fourier system. Should you wish to alter the allocations, you will also have to 
modify the keyboard programs to reflect such changes. 
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POWER SPECTRUM Wadena nos 
FLOW DIAGRAM (STACK 57) 


SET: AOC SWITCHES - 
BLOCK SIZE 


(1.1100) 


(L1101) 
3 SET HP FILTERS 


LABEL 2000 ; 
SET # AVERAGE 









(11250) 
(L1252) (L1251) 
SELECT TYPE BSFA 
AND CF/BW ne 
(1300) BASEBAND 
NUMERICAL 
x ga 
(11301) (11302) 
ENTER CF/BW SET BW 
FROM TERMINAL WITH CURSOR 
(3000) 
(1450) 
BASEBAND TRANSFER 
OLO Te NEW FUNCTION 


MEASUREMENT 





(L1553) 


SET THRU—PUT 
TRACKS 
a (1600) 


ON—LINE PRE—PROCESSOR 
OFF-LINE pel 
SOFTWARE TYPE 
, - BSFA 
? OFF-LINE 


PRE—PROCESSOR 
(L1603) (L1601) 


INITIALIZE BSFA INITIALIZE BSFA INITIALIZE BSFA 








OFF-LINE OFF—LINE ON-LINE 
SOFTWARE BSFA 54470A BSFA 54470A BSFA 







(L1259, L1609) 


an DISPLAY BLOCKS & 
PLOT IF DESIRED 
PRESS CONTINUE Gori > 






— 


—_—- ee 


TURN ON 
RANDOM SOURCE 






TRANSFER FUNCTION 
FLOW DIAGRAM 


L1000 
LABEL 1000 
(STACK 51) 
INPUT TYPE OF 
TRANSFER FUNCTION 


(L1050) 
(L1052) 
YES 
INPUT DAC VOLT. 
(L1051) 
(L1053) 
NO 
(L1100) 
INPUT BLOCK SIZE & 
FILTERS (IF AVAILABLE) 
(L1150) 


c 


TURN ON HP DAC 






(L1153) 


SET TRIGGER SOURCE 


SET UP ADC SWITCHES 
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LABEL 2000 
SET # AVERAGE 


(L1200) 


TYPE 
TRANSFER 
ENS 


(L1202) 
TURN ON HP DAC 










RANDOM X 











BSFAN < mer BASEBAND 


SELECT TYPE BSFA 
AND CF/BW (L3000) 


(L1300) BASEBAND TRANSFER 


NUMERICAL FUNCTION 


CURSOR MEASUREMENT 


ENTER CF/BW SET BW YES WAS 
FROM TERMINAL WITH CURSOR DAC ON 
? 





OLD TYPE NEW 






(1500) 


(L1553) 






SET THRU—PUT 
TRACKS 


USE LISTEO 
TRACKS 


(L1600) 


ON-LINE PRE—PROCESSOR 
OFF-LINE rh 
SOFTWARE TYPE 
BSFA 
f OFF-LINE 











PRE—PROCESSOR (L1601) 


INITIALIZE BSFA 
ON-LINE 
54470A BSFA 


(L1603) (L1602) 


INITIALIZE BSFA INITIALIZE BSFA 
OF F—LINE OFF~—LINE 
54470A BSFA 





SOFTWARE BSFA 


(L1209, L1259, L1609) 


DISPLAY BLOCKS & 
PLOT IF DESIRED 


PRESS CONTINUE <> 
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END OF BLOCK £100 
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BLOCK £150 
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END OF BLOCK 14150 
Analog input in "REPEAT" to monitor inouts 


Power spectrum (log) for user reference 
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STACK 53 


Cursos on, Cursor parameters to 


variables 


Variable 


Cursor on, parameters to 2000 - 


Cursor 


0 
Variable 


2000 - 2002 (200L=frequency) 
2003 gets ist cursor frequency 


2002 


ff 
2004 gets 2nd cursor frequency 


BLOCK 1356 
If ist cursor frequency ) end 
then swap thes 


END OF BLOCK 135) 


Variable 


20 gets zoom bandwidth 


Variable 19 gets 
BLOCK 1400 
If galcucdet not { ctr freq 


then bandwidth = ctr freq - 


END OF BLOCK 140i 
END SF BLOCK 1200 | 
BLOCK 1450 (ends in stack 54) 

Reranch per yvariadle :7 

to L£4Si, 1452 or 1453 


iTPom stack 52) 


Zoom center frequency 


2004) 
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Variable: 24 
BLOCK 1500 (ends in stack 54) 
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STACK 54 


Var.able 22 gets default start track 
Variable 2f gets default ¢ of records 


Parcmeter entr, 
Range: {£ - 197 
Variable: 22 
Variabie 22 qets yser’s start track 
yoriable Ci gets ¢ af recerds iert 
END CF BLOCK 1550 .from stack 3? 
Position throvgnput file te start track 


Pouse for snerator act.on 

Variagie £3 gets current Glock 3ize 
Set olock Size 79 aay for tarougnput 
Perform ADC throughps% 

Restore clack size 


END GF BLACK £599 (frsa stacx 53) 
END OF BLUCK 2459 (fron stack 53) 
BLOCK £480 


Branch per yariable i? — 
to Liog{, 1602 or idds 


Thizagiize roomMm~on L228 oreprocessar 
Subroutine call--zeom measurement 
Reset room 76 Daceband 

Ge to end of block isit 


initialize zoom--etf line prepracessar 
Subroutine call--zoom# seasuresent 
Reset zoom te basebanc 

Go to end of biock idté 


roame-3o4? Line software 
i Subroutine call--zo0m measurement 
Reset zoom te dase2and 
END OF BLOCK iSd0 sfrim stack 49: 
END JF SLOCK i257 ifrsm stacx $2) 
END OF BLOCK {229 (from stack $2) 


Transfer function cimplete 
Print .nstructions 10 operator 


. 
“Nitzalize 


Save variable paraneters 


Pause for operator action 
_, LOD to eyed transfer function 
END OF TRANSFER FUNCTION (except suaroutines) 


en 
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“55 STACK SS : 

as BASERAND TRANSFER FUNCTION SUBROUTINE 
6 2 3 Check ADC input selector. If not 2 channel 
3 4 inform dperator 


Quse for operator to correct it and 
3008 =H oop to check again. 


: Clear blocks needed for averacing 
4 e 
5 
5050 BLOCK 3250 - 
J 3050 {1D Branch oer variable Li 
_. 0D ={ te SBS, duce cp 305s 
ate L305 & 3052 equivalent--nan-transient 
3 Analog input 
q Hann 
{ 
i Transfora 
9 j Clear DC o 
r Q 
d 2 : Averace spectrs 
sect iS5 ) Laop ror specified # of qvernaes | 
7659 7 1 Subroutine cali - clear 2nd nal? diocs 
t 2 é Compute transfer function - 
fe to Correct ampiitude for hanning 
% 3 
3 {5 . 
3 3 
3059 =i . 
3053 L3053 transient siancls 
2 { Instruct operator ta capact on cue 
3054 . 
.) wue with “bees” 
5 Analog .apet ; 
S J =u Check for sverlcad--it so, cue again 
q i Transfore. : 
q J Ciear iC 
$ Q 
Q e 2 Sverace spectra 
3554 1S) 0 Loco cor cdecaficd # of averages 
9059 fi i Subesutine cai: - cigar 2nc half sloc: 
0 2 Z ‘Cemoute thansrer runction 
3059 END CF BLOCK 3650 
: yake 1098 of spectra 
3 
4 


Retura to calling orecten 
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fo ae | FO 


STACK So 
SUBROUTINe SOR ON LINE ZOUM TRANSFER FCN 


Check ADC input selector If not 2 channels 
inform operator 
pause for operator <9 correct it and 
loop to check again 

Clear dlocks needed for averaging 


Zoom, displaying Channel & power Spectr 
Averace spectra 

Loop for specified # 2f averages 
Comdute transfer function 

Subroutine cail te move cata 

Return tc calling pregran 


SEP SS SS SF SBF SOF S88 S82 228 8 SB S88 “SB 2 288 S282 @2e 8 SS = = 
EHS Seb A p> ee ep ee rm ae eG mm ee 


SUBROUTINE 10 A0Ve DATA BLULKS 


Move data DiacksS from room area 
to area for operator viewing 


Take Logs of spectra 


Return to calling aregram 


Set OES OS SEES SES SESE SS SSP SS SCE SP 2882 Se 2 eS eS S we 
abel a> a> GD ae ae ab US OD oD Oe OD ee ee es ee 


SUBROUTINE FOR OFF LINE ZOOM TRANSFER rCN 


Clecr olocks needed for averaqing 


Zoom on Channel 
Zoom on Channel 
Display Channel B sower spectrum 
Average spectra 

Loop ror specified # of averages 
Compute transfer function 
Subroutine call to nove data 
Return to calling progres 


a> rad 


ee ee a oe a Oe 2 SS BSS SS SSS SSE SSB SF S88 BSB SES SS OSs 
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STALK $? 
POWER SPECTRUN PROGRAS 


Recali Previous vericole puraseters 


Read text buffer #Si 
Print messages 445 


Blacn size Si2--yser aqy change after pause 
Pouse for operator action 


Parameter entry 
Rang]: « 7 4 
Varzabie: 34 


BLOCK £209 


307 2 SLbtersm-quto selert méae 


END OF RLOCK 4496 


Peuse for operator actin 


mhclog input in "REPEAT" t9 monitor snouts 


Power spectrys 140g) rom uSer rererenc: 


Jump to next stack and continue 
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FF 


STACK 58 
<3 LOOP POINT FOR POWER SPECTRUM 
{ 


{ Parameter entry 
52767 Range: 4 - 32767 


0 Varsable: 15 
atti 
{ Parameter entry 
z Range: 1 - 2 
0 Variable: 16 : 
2000D 


BLOCK £250 (ends in stack 50) 


{250 16D Branch per variable 46 
VOmeteoiedn tode 


Pouse for operator action 


3 { Subroutine call - baseband measurement 

2 { then go to end of block i 
{ 
{ Parameter entry 
3 Range: i - 3 

0 Variable: 17 

2000D 

aS 
{ 
{ Parameter entry 
2 Range: { - 2 

0 Variable: 48 

2000D 


BLOCK 4300 (ends in stack 59) 
{300 i8D Branch per variable 48 
to Lig04 or 4302 


Parameter entry 
32767 Range: { - 32767 


0 Variable: i7 
2000D : 
i . 
ie | Parameter entry 
17) i Range: 1 - <value of var 19 - 4) 
0 Variable: 2 
26000 
i { 
i i 
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STACK 59 


Cursor on. Cursor parameters to 
variables 2000 - 2002 (2001=frequency) 
Variable 2003 gets ist cursor frequency 


Cursor on, parameters to 2000 - 2002 
Cursor off 
Variable 2004 gets 2nd cursor frequency 
BLOCK 1350 
If {st cursor frequency ) 2nd 
then aa thes 
END OF BLOCK 1350 


Variable 20 gets zoom bandwidth 


Variable 19 gets zoom center frequency 
BLOCK {400 
If dandwidth not ( ctr freq 
then bandwidth = ctr freq - { 
END OF BLOCK 1400 
END OF BLOCK 1300 .fram stack 58) 
BLOCK 1450 (ends in stack 60) 
Branch per variable i7__ 
to Li4Si, 1452 or 1455 


17D 
Pause for operator action 


Parameter entry 
Range: { - 2 
Variable: 24 

BLOCK £500 (ends in stack &0) 


Branch per variable 24 
to LiSOL or 1502 


24) 


Parameter entry 
ange: i -_ 3 
Varzable: 2000 (temporary) 
BLOCK 1550 (engs in stack 605 
Branch per variable 2008 
te LIGSi, i552 or 1553 


"20000 


Pause & loop te here--abort 
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2000D 


{7D 


STACK 60 


Variable 22 cets default start track 
variable 2{ gets default € of records. 


Paraneter ante. 
Range: 1 - {77 
Variable: 22 

Variable 22 gets yser’s start track 
Variable 2f gets # of records left 
END QF BLOCK 1550 (from stack 59) 
Position throughput file to start track 


Pause for eperator action 

Variable 13 gets current block size 
Set block size to max for throughput 
Perform ADC throughput 

Restore block size 


END OF BLOCK {500 {from stack 59) 
END OF ELOCX 41450 (from stack $9) 
BLOCK 1660 


Branch ger variable {7 
to Li604, 1602 or 1463 


Initialize zoom--on line preprocessor 
Subroutine call--zoom measurement 
Reset zoom t4 Daseband 

Go to end of block 1604 


Initialize zoom~-off line preorscessor 
Subroutine call--zoom measurement 
Reset zoom to baseband 

Go to end of block 1604 


Initialize zoom--uff line software 
Subroutine cali--zo0m seasurenent 
Reset zoom to baseband 
END OF BLOCK 1600 
END OF BLOCK i250 (from stack $9) 


Power spectrum complete 
Print messages to operator 


Save variable paraneters 
Pause for operator action 


ie to repeat ae spectrua 
END OF POWER SPECTRUM (except subroutines) 
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STACK 21 
BASEBAND POWER SPECTRUM SUBROUTINE 


Check ADC input selector. [f aot 1 channel 
i inform operator 
Quse for operator to correct it and 


fo§ 


=i oop to check again 
Clear block needed for averaging 
Anaiog input, displaying average 
ann 
Transfors 
0 Clear DC 
Average spectrun 
{52 0 Loop for specified # of averages 
t { Subroutine cali - clear end half black 
1§ Hanning correction for broadband noise 
Load average to block J 
Take log 
Return to calling program 
SS “SUBROUTINE FOR ON LINE ZOUH POWER SPECTRUM 
{ 3 
i Check ADC selector. If not 2 channels 
. inform operator 
a pause for operator to correct it and 
.. loop to check again . 
Llear biock needed ror averaging 
{ | 


Zoom, displaying power spectrum 
15D 0 Average spectrum 

Loop ¢or specified # af averages 
Load average to block J 
Take log 


ttn et ast eS Ss KM SS SO SOB SES BPs Fe SB Sees ee Ss ee 
smb cmb ont ans GD aes GD OES GD GE aes es ee ee ee ee ee ee ee ee eae eas a a ee 


SUBROUTINE FOR OFF LINE OOM POWER SPECTRUM 


if Clear dlock needed for averaging 
4 
Zoor 
19D 0 Display averace ({ sweep) 
Average 


Loop tor specified # of averages 
Load average to dlock J 

Take Lo 

Retura fo calling orcegran 
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L “a2 STACK 62 

L 9050 CLEAR 2ND HALF BLOCK SUBROUTINE 
ie 27 057 Set return label valve 

J 3056 -{ 

L 9057 Return label 

CL & 25D 26D Clear last half block 4 
CL S 25D 26D Clear last half block $ 

J 9059 ={ 

L 9055 Entry pt. Pwr Spec a 
em 2? =69058 Set return label valve 

J 9056 -{ 

L 9058 Return label 

me { 25) 26D Clear last half block 4 

J 9059 -{ 

L 9056 . 

Y at Q 905i {4D Wes dP Filter selected? 

J 0D -{ No, Gto 9051 

Eee 7 052" Yes 

Y BS i Get current blocksize 

ie fd {3} Store Bs/2 in variable 26 
ee 2 26D Store 8s/4 in variable 25 
J 29) =y Indirect return 

tL = s- 054 

L 9059 End of ogm 905049055 

( Return to calling progran 


Text buffer messages for Transfer Function and Power Spectrum programs 


BUFFER # MESSAGE # 
Si { 
HP TRANSFER FUNCTION PROGRAH 


Sf 2 
SELECT EXCITATION TYPE 
S{=RANDOM - RASERAND OR ZOOM 
2=HP DAC - BASEBAND ONLY 
S=TRANSIENT - BASEBAND ONLY 


St g 

INPUT DESIRED DAC GUTPUT IN MV 
5f 3 

HP POWER SPECTRUA PROGRA 


St 5 
SET ADC FREQUENCY RANGE AS DESIRED 
(SAMPLE MODe & MULTIPLIER) 
SET ADC TRIGGER (0 “FREE RUN" 
CHANGE BLOOK SIZE iF DESIRED 
PRESS "CONTINUE" WHEN READY 


51 
ARE HP FILTERS INSTALLED? 
O=H0, Si=YES 
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Si 7 
oy sa REPEAT/SINGLE SWITCH 
PRESS “CONTINUE” WHEN READY 


ot 98 
UPPER LIMIT = 


St ve 
ee ee RANGE-PLEASE REENTER 
PRESS CONTINUE WHEN EADY 


“ye 8 
TURN ON RANDOM EXETTATION SGURCE 


Se 7 
SET TRIGGER SOURCE AS DeSIRED 
IMPACT STRUCTURE REPEATEDLY 


Se 10 
SET JVERLOAD VOLTAGES AND TRIGGER 
LEVELS FOR SIGNAL ar ATU 
AOVe REPEAT/SINGLE SWITCH 72 
TSTNGLE” WHEN READY. iF SQURCE 
TS NOT IN FREE RUN. TRIGGER THE 
SYSTEM AGAIH 10 CORTINUE. 


Be tf 
ENTER NUMBER OF AVERAGES DESIRED 


$2 ce 
ENTER MEASUREMENT 7 
{=BASE RAND 
52=200d 


YPE 


Se {3 
ENTER Z00M MEASURENEAT MODE 
{=0N LINE, PREPROCESSOR 
S2=OFF LINE, PREPROCESSIR 
3=0FF LINE, SOFTWARE 


149 





Se 14 
ZOOM NOT APPROPRIATE WITH HP DAC 
BASEBAND MEASUREMENT WILL BE MADE 


$3 ue 
HOW WILL YOU SPECIFY ZCO0K ENE 
{=NUMERIC ENTRY - CTR FREQ & 
2=CURSOR ~ ON PRIOR MEASURERENT 


53 {6 
ENTER CENTER FREQUENCY 

53 7, 
ENTER BANDWIDTH 


53 13 
MOVE CURSOR TQ START FREQUENCY 
PRESS "VALUE (SWITCH REGISTER {13 


33 1D 
MOVE CURSOR TO END FREQUENCY 
PRESS *VALUE" 


$3 20 
ANALYZE OLD OR NEW DATA? 
ear (FROM THROUGHPUT FILE) 


54 ag) 
THROUGHPUT COMPLETED 
5 oe 
PRESS "CONTINUE" FOR MEASUREMENT 


54 a3 
MEASUREMENT COMPLETE 


C4 
TO DISPLAY RESTS PRESS: 
"DISPLAY® "0° LGG TRANSFER FCN 
"DISPLAY" "4" COHERENCE 
"DISPLAY" *2" INPUT POWER SPECT 


50 


Beeps to cve operator 





"DISPLAY® "3" QUTPUT POWER SPECT 
"DISPLAY" °54" CROSS POWER SPEC’ 


2 


54 
TO COPY DISPLAY. ON TERMINAL: 
PUT TERMINAL IN GRAPHICS NUDE 
PRESS “GOLD KEY" "PLOT" 


54 26 
TG MAKE ANITHER MEASIREMENT : 
PUT VERMINAL IN ASCE. MODE 
PRESS “CONTINUE 


2) 
= 
| 


54 7 
ENTER STARTING TRACK FOR THRGUGHPUT 


G4 ce 
THROUGHPUT WILL USE TRACKS 
£35 THRGUGH 199 GN THE LOWER 
(FSDS) DISC. IS THIS UK? 
‘SG NO - ABORT 
e=YeS - PROCEED 
3 NO - ASK ME FOR TRACK $ 


55 ay 
THPACT STRUCTURE ON CUE <REEP) 
FOR EACH AVERAGE 


on 30 


55 at 
SET ADC INPUT SELECTOR TO "AR" 
PRESS °CONTINUE* WHEN Ready 


5S 32 
SET ADC INPUT SELECTOR TO *A" 
PRESS “CONTINUE* WHEN READY 


eel 


Reens to 


Beeps to 


Beeps t0 


cue operator 


cue operator 


cue operator 
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APPENDIX C 
CALIBRATION DATA 


acl ak Se) ee 


E ° 


t 


ACCELEROMETER MOUNTING WAX 
Model OQ80A24 





PIEZOTRONICS 


Petro-Wax functions to transfer motion from the test object to the sensor. It is 
used to couple the sensor directly to the test object. It is a convenient, temporary mount- 
ing replacement for studs and permanent adhesives. 

An elastic interface structure with heavy damping, it forms quickly to irregular sur- 
faces to facilitate easy mounting of sensors. This pliable wax enables the sensor to be 
mounted in nearly any convenient spacial coordinate system. 

In conjunction vith tine sensor, the wax Forms a heavily damped spring-mass system, 
The frequency response is a function of transducer mass, mounting area, depth of wax, and 
test temperature. As the amount of wax at the interface of the sensur and test object in- 
creases, the first resonance of the system decreases, which limits the frequency response 
of the fixturing (see graphs). - This mounting technique is primarily for use at room temp- 
erature. The wax cannot be used effectively at high or low temperatures. 


How to use: (1) To insure a secure and valid fixture, be certain all surfaces are Free 
of oll and dirt. 

(2) Apply the wax directly to the base of the transducer or to the adhesive 
mounting base. The amount of wax used depends on your individual appli- 
cation. Use of the adhesive mounting base helps keep the transducer clean. 

(3) Press the sensor firmly against the test structure to insure secure mount- 
ing and as little vax at interface as possible. [his provides the best 
frequency response. . 

(4) Procede with measurement. 

PETRO-WAX is availuble in quantity from: Katt & Associates 
E.G. Hox. I8269 
Pittsburgh, PA 15227 
(412) 885-5727 








| 8B 

MODEL NO >, Thine | <0;010" 

Horizontal Scale % UOO Hz Magnitude 

Mag. : vert. scalel0 db/div Phase | 
Phase: vert. scale 50 °/div Frequency a 





Thin Med Thick 





Med “0.030" 
Magnitude 


Thick 950.125" 
Magnitude 


| 500 | 10 | 
fio; || 
20 000 ra 


PCB PIEZOTROMICS, INC. 34e5 WALDEN CVENUE DEPEW. NEW YORK 14043-2895 TELEPHONE 716-684-0001 TWX 710-263-1377 





Ue 3) 





MODEL NO. 302A 
Horizontal Scale 10 000 Hz 
Mag. : vert. scale 10 db/div 


Phase: vert. scale 50 °/div 


Thin Med Thick 





MODEL NO. 303A 
Horizontal Scale 25 000 Hz 
Mag. : vert. scale 10 db/div 


Phase: vert. scale 20 °/div 


Thin Ned Thick 





Thin <o.010" 


Magnitude 


Frequency —> 
t 











Med ~ 0.030" 


Magnitude 


Phase | ioe 
| oe lo KHz 
Thick >0.125" 
Magnitude 
Phase 
_ 10 KHz = 


Thin <o.010" 


Magni tude 





! 


| Med “J .U30" 


Magnitude 


Thick pa OP bee 1 


Nagnitude 


1 





PLB CALIBRATION CERTIFICATE - BER 
MOACT FORCE HAMMER i ui 
ho. : ws FIEZOTRCNICS, INC. 
Model No. C8568 O38 ~. 

7 = Customer: 
Transducer Model No. — ; 2 E 

° “ : Ad vas. STeqUALIATA PC owas 
Transducer Serial No. 2067 Ss a 
Hammer Calibrator: PCB Model Zea Invoice No: 
Pendulous Mass_ 2, 7-1 (JP gm) oe 


including accelerometer 


: 
—_—a=—— 


Initials: fe 1 


Basic Force Transducer Sensitivity JO bE mV/1b (reference) 


(for stationary installations) 


Cate: mae 3 mas 


Reference Accelerometer Sensitivity JE, | - mV/Q (302A'%/N Fo) 

-— ; a \ 
Ratio: force Transducer Sensitivity/Accelerometer Sensitivity {OS | (Gage 
(from test of mass impacting stationary hanmer) 


HAMMER SENSITIVITY (3) (with 2.7 gm Al/Plastic, Al/Steel, Al/Rubber 7:95; 


with with with 
Steel Extender | Al. Extencer | Wo Extencar 


we ae 
Cc GY 7 oe 


\ 
ase | [Pee 
Pere aNE 


HOTcS: (1) Difference from reference sensitivity of force transducer 


Configuration 









_ —<—- ~-= 






Ratio: Force/Acceleration 
Sensitivity (2) 


-—-—_— — a | 


Hammer Sensitivity mV/1b 


Difference: (1) 


(2) In transfer fuction testing, the important factor is the ratio of 
sensitivities. 


(3) Because of normal behavior of the hammer structure, the aoparent 
sensitivity of the hammer in motion differs from the stationary 
calibration of the force transducer. [It is less by a factor 
proportional to the ratio of the ness of the imoact cap and Seismic 


se distributing mass in the transducer to the tctal mass of the 


henvaer structure. Using a heavier hammer head or installing the 
force transducer on the structure cnly chanaes the problem. A 
—--- heavier hammer head tends to penetrate the test object or cause 
multiple bouncing. «hen mounted on tne test object, the inertial 
mass in the transducer causes it to act as an accelerometer sensing 
motion of the test object. 


ro4 









Bes tiv, Vere eho ee Model No. _ 32 3.79 S 





| CALIBRATIU:. .TA 
sed * for Serial No. _ 13 4 ; 
i 1.C.P.ACCELEROMETER | Range $d 





(per 1S A $37.2) 
PIEZOTRONICS 


Max Input oe 4 i 


ie ND aoe Max Teinp fe 


Meee. ESM @ 100 Hz,8q9sPK 


2. MAXIMUM TRANSVERSE SENSITIVITY Me Z percent 


Order No. NGZ271183 © oudeE 
ana JA bo 





1. VOLTAGE SENSITIVITY 











6. 
toate | | to | 20 | so | 00 | 200 | soo | somo | soso | sro | Zree pave 
[owiiont | 23h 7.0 | 


Calibration traceable to NBS through project no. 






PCB PIEZOTRONICS. INC. 
3425 WALDEN AVENUE 
DEPEW. NEW YORK 14043 


Customer Mae Veter cware Gee Model No. x OA LA a Z 


CALIBRATION DATA 


tae SerialNo. 2 P7 
; 1.C.P. ACCELEROMETER Range ec ¢ 
; , (per tS A $37.2) 


Max Input oe aId ‘ 
ony, }¢ I! Max Temp Cc Se C 


1. VOLTAGE SENSITIVITY ...... f4. [A mv/g }| @ 100 Hz,8gsPK 


Order No. .74, 2-277) B34 F ove, 





2. MAXIMUM TRANSVERSE SENSITIVITY EG re 
3. RESONANT FREQUENCY ......... ce fe Sete 
4. DISCHARGE TIME CONSTANT ..... O.F peeonts 


5. OUTPUT BIAS LEVEL 


6. 


Pome |b ebeg baal + a0 basen 


ae 


3425 WALDEN AVENUE 


3 
OEPEW. NEW YORK 14043 ay = ne 
3 v7. 


oe ee @ ee @ @ @ @ 


FREQUENCY RESPONSE: 
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IMPULSE FORCE HAMMER. 


Hodel No. 0862 FO Kit No, GRAY SE Cust omer: " 
Serial No. 2s5 : hw Ee he fe Airy wa TH ee ite iw 
Range F- S% kb | 


Linearity error ~A.e & Invoice No,: 


Discharge Time Constant_ “7° 8 ALU A 
Output Impedance /C¢ ohms 

Output Biss /Z.2, volts 

Traceable to N8S through Ves dase 

Initials “ACKL pete: 4 /A/ 83 


FoRce 
Aocelerameter: Mode] No. ZO38 __, Serial no. 2974 , sens. C72 mV/q 


Pendulous Test Mess —— 1p cram) including scceleropeter 





Hammer Sensitivity: 


CONFIGURATION 


Extender 


HAMMER 


SENSITIVITY 





NOTES: 


1. The sensitivity ratio (Sa/Sf) is the scsling fsctor for converting structural transfer 
messurements into angineering units. Divide results by this ratio. 


2. Each specific hammer configuration hes a different sensitivity. The difference is 


@ constant percentage, which depends on the mass of the cap and tip ascendly relstive 


to the tots] macs of the head. Calibrating the specific hammer structure being used auto 
matically canpensates for mass effects, 


Effective mass [— |} with 302ACQ7 attached and vinyl capped plastic tip. 
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eK I.C.P.ACCELEROMETER | Range GO? 1: 
:.,, ee : PIEZOTRONICS . eae Max Input__§ “7°77 a 
) . 150 bau Nemantec Me Max Temp __ BY OF 





1. VOLTAGE SENSITIVITY ...... S, bs mv/g @ 100 Hz,8gsPK 


2. MAXIMUM TRANSVERSE SENSITIVITY a, O ” percent 


3. RESONANT FREQUENCY ......... eee a 


4. DISCHARGE TIME CONSTANT ..... Aart seconds 


5. OUTPUT BIASLEVEL ........... _ &e 7 Volts 


PROCUVENUT Neoruivoc. 


ee NS sd 
Lowi | EERE AGE / | a [paplecalecPless |p3.ale 3 |50) 


Calibration traceable to NBS through project no. Bip isaae 


PCB PIEZOTRONICS. INC. ess? a 


3425 WALDEN AVENUE : 
DEPEW. NEW YORK 14043 a 


Model No. OS A7OT_ 
CALIBRATION DATA 
for Serial {e927 
1.C.P. ACCELEROMETER 
(per 1S A $37.2) 
PIEZOTRONICS 





















tomer VA = 
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Max Temp ae OF 


5t + -th le bot 


UV VOLTAGE SENSITIVITY ....4.. 1/: PO mug @ 100 Hz,8 gs PK ’ 


2. MAXIMUM TRANSVERSE SENSITIVITY / 7, percent 
3. RESONANT FREQUENCY ......... AGS KHz 


. FREQUENCY RESPONSE: 


[ele [al clolelelebeles 
Powe | be alystral » (vaAciolnceealerdszdez/ 


Calibration traceable to NBS through project no. 42 Eee 


PCB PIEZOTRONICS, INC. wt / et 


3425 WALDEN AVENUE 
DEPEW, NEW YORK 14043 
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APPENDIX D 


LOCALLY GENERATED USER KEYBOARD PROGRAMS 


The locally generated keyboard programs presented in the 
next six pages allow the use of impulse technique or random 
excitation while utilizing the Modal package of the HP-5451C 
Fourier Analyzer. Each program and subprogram is identified 
by the first line of the program listing. Program number 1 
(identified by the -1l on line 1) is the control program, 
feet rom this program subprograms 50, 51, 52, 53, 54, 58 
and 59 are called as needed. 

All eight programs are automatically loaded into the 
computer memory when the following command is executed, as 


discussed in Chapter V, step number 35, 


[JUMP] 0 1 [ENTER] 
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IF 


2 -2 
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-2 
2 0 
2 0 
9 0 
1 0 
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125 
iol 
135 
139 
145 
151 
155 
1635 
169 
175 
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189 
195 
199 
206 
eile 
217 
224 
250 
237 
244 
251 
258 
262 
269 
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280 
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292 
2o7 
303 
311 
315 
321 
324 
328 
5354 
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346 
351 
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290 2 
290 30D 0 
2 0 2 
i 0 i 
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210 0 
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1 
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0 2 5606 
1 
0 
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2 2 
42D 410 
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2 rs 
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MODE CHz) €X%) CRAD/SEC) 
emerges =6hlllCtC<C~mt”St‘(Cti‘sS™S™S™S™S™*~:”~<~iR 
| pt gBd3.allt! 0.4564 i: One” 

| | t 
Z| 2362-30} 1.6598 | 246.3933 | 
| | | | | 
me) gecoego! Oss ! 128.9980 | 
1 | t | 
| Y | 24.7603 | 9.q5qo | 69.2302 | 
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